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Abstract 
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1 Introduction 

There is a number of good textbooks in quantum field theory (QFT, see [U El [3l HI [U El [71 
[HI [HI [To]. However, they are all designed for the graduate- level study and we can only hope 
that undergraduate students can read some easy parts of them. Moreover, there are certainly 
no undergraduate-level textbooks for TQFT, so pure students are forced to try to read the 
original papers from its inventors, Ed Witten [111 I12j and Michael Atiyah [iBJ. The goal of 
the present lecture notes is to try to fill in this gap, to give the talented undergraduates the 
very first glimpse of the mathematical physics of the XXI Century. 
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Throughout these lecture notes we will use the following conventions: 
(i) natural units, in which (some or all of) the following definitions are used: c = h = m = 1; 
{u)i = y/^, z = dz/dt, dz = d/dz; 

(iii) Einstein's summation convention over repeated indices, while nD means n— dimensional. 
1.1 Basics of non-relativistic quantum mechanics 

Recall that quantum theory was born with Max Planck's 1900 paper, in which he derived the 
correct shape of the black-body spectrum which now bears his name, eliminating the ultravi- 
olet catastrophe - with the price of introducing a 'bizarre assumption' (today called Planck's 
quantum hypothesis) that energy was only emitted in certain finite chunks, or 'quanta'. In 
1905, Albert Einstein took this bold idea one step further. Assuming that radiation could 
only transport energy in such chunks, 'photons', he was able to explain the so-called photo- 
electric effect. In 1913, Niels Bohr made a new breakthrough by postulating that the amount 
of angular momentum in an atom was quantized, so that the electrons were confined to a 
discrete set of orbits, each with a definite energy. If the electron jumped from one orbit to a 
lower one, the energy difference was sent off in the form of a photon. If the electron was in 
the innermost allowed orbit, there were no orbits with less energy to jump to, so the atom 
was stable. In addition, Bohr's theory successfully explained a slew of spectral lines that had 
been measured for Hydrogen. The famous wave-particle duality of matter was proposed by 
French prince Louis de Broglie in 1923 in his Ph.D. thesis: that electrons and other particles 
acted like standing waves. Such waves, like vibrations of a guitar string, can only occur with 
certain discrete (quantized) frequencies0 In November 1925, Erwin Schrodinger gave a semi- 
nar on de Broglie's work in Zurich. When he was finished, P. Debye said in effect, "You speak 
about waves. But where is the wave equation?" Schrodinger went on to produce and publish 
his famous wave equation, the master key for so much of modern physics. An equivalent 
formulation involving infinite matrices was provided by Werner Heisenberg, Max Born and 
Pasquale Jordan around the same time. With this new powerful mathematical underpinning, 
quantum theory made explosive progress. Within a few years, a host of hitherto unexplained 
measurements had been successfully explained, including spectra of more complicated atoms 
and various numbers describing properties of chemical reactions. For more details, see, e.g. 



1.1.1 Quantum states and operators 

Non-relativistic quantum-mechanical systems have two modes of evolution in time [T6] . 
The first, governed by standard, time-dependent Schrodinger equation: 



^The idea was so new that the examining committee went outside its circle for advice on the acceptability 
of the thesis. Einstein gave a favorable opinion and the thesis was accepted. 
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describes the time evolution of quantum systems when they are undisturbed by measure- 
ments. 'Measurements' are defined as interactions of the quantum system with its classical 
environment. As long as the system is sufficiently isolated from the environment, it follows 
Schrodinger equation. If an interaction with the environment takes place, i.e., a measurement 
is performed, the system abruptly decoheres i.e., collapses or reduces to one of its classically 
allowed states. 

A time- dependent state of a quantum system is determined by a normalized, complex, 
wave psi-function ip = il'{t). In Dirac's words tl5j, this is a unit ket vector which is an 
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element of the Hilbert space L'^{tl^) = Ti, with a coordinate basis The state ket-vector 

\ip{t)) is subject to action of the Hermitian operators, obtained by the procedure of quanti- 
zation of classical biodynamic quantities, and whose real eigenvalues are being measured. 

^The family of all possible states 10), etc.) of a quantum system confiture what is known as a Hilbert 
space. It is a complex vector space, which means that can perform the complex-number-weighted combinations 
that we considered before for quantum states. If lip) and j(j!)) are both elements of the Hilbert space, then 
so also is w\il}) + z\(j)), for any pair of complex numbers w and z. Here, we even alow w = z = 0, to give 
the element of the Hilbert space, which does not represent a possible physical state. We have the normal 
algebraic rules for a vector space: 

W + \<l,)^\ct'} + W, 

M + m + \x))-m + m + \x), 

w{z\%p)) — {■wz)\%p), 
(to + z)\rp) = wlip) + z\tP), 
ziW + m ^ z\i<) + z\^) 
0|V') = 0, zO = 0. 

A Hilbert space can sometimes have a finite number of dimensions, as in the case of the spin states of a 
particle. For spin |, the Hilbert space is just 2D, its elements being the complex linear combinations of the 
two states j t) and | J,). For spin in, the Hilbert space is (n + 1)D. However, sometimes the Hilbert space can 
have an infinite number of dimensions, as e.g., the states of position or momentum of a particle. Here, each 
alternative position (or momentum) that the particle might have counts as providing a separate dimension 
for the Hilbert space. The general state describing the quantum location (or momentum) of the particle is 
a complex-number superposition of all these different individual positions (or momenta), which is the wave 
'i/'— function for the particle. 

Another property of the Hilbert space, crucial for quantum mechanics, is the Hermitian inner (scalar) 
product, which can be applied to any pair of Hilbert-space vectors to produce a single complex number. To 
understand how important the Hermitian inner product is for quantum mechanics, recall that the Dirac's 'bra- 
kef notation is formulated on the its basis. If we have the two quantum states (i.e., Hilbert-space vectors) are 
\^) and then their Hermitian scalar product is denoted {tp\4'), and it satisfies a number of simple algebraic 
properties: 

{'4'\4'} = {(f'Wj (bar denotes complex-conjugate) 

{zi-^m^zim, 

W^) > 0, i^pw = if i^) = 0. 

For example, probability of finding a quantum particle at a given location is a squared length \tp\'^ of a Hilbert- 
space position vector which is the scalar product of the vector |i/)) with itself. A normalized state 
is given by a Hilbert-space vector whose squared length is unity. 

The second important thing that the Hermitian scalar product gives us is the notion of orthogonality between 
Hilbert-space vectors, which occurs when the scalar product of the two vectors is zero. In ordinary terms, 
orthogonal states are things that are independent of one another. The importance of this concept for quantum 
physics is that the different alternative outcomes of any measurement are always orthogonal to each other. 
For example, states | t) and | I) are mutually orthogonal. Also, orthogonal are all different possible positions 
that a quantum particle might be located in. 
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Quantum superposition is a generalization of the algebraic principle of linear combination 

of vectors. The Hilbert space has a set of states \^p^) (where the index i runs over the degrees— 
of-freedom of the system) that form a basis and the most general state of such a system can 
be written as 1^^)= ^i^i Iv'j)- The system is said to be in a state \i^{t)), describing the 
motion of the de Broglie waves, which is a linear superposition of the basis states \(fj} with 
weighting coefficients q that can in general be complex. At the microscopic or quantum 
level, the state of the system is described by the wave function \tp) , which in general appears 
as a linear superposition of all basis states. This can be interpreted as the system being in 
all these states at once. The coefficients q are called the probability amplitudes and |cj| 
gives the probability that \tp) will collapse into state \ip) when it decoheres (interacts with 
the environment). By simple normalization we have the constraint that |cj|^ = 1. This 
emphasizes the fact that the wavefunction describes a real, physical system, which must be 
in one of its allowable classical states and therefore by summing over all the possibilities, 
weighted by their corresponding probabilities, one must get unity. In other words, we have 
the normalization condition for the psi-function, determining the unit length of the state 
ket-vector 

imim) = j r^dv = j ivi'rfF = i, 

where V'* = (^(01 denotes the hra vector, the complex-conjugate to the kct ■(/' = ^-^(0)) ^iid 
{il){t)\il){t)) is their scalar product, i.e., Dirac bracket. For this reason the scene of quantum 
mechanics is the functional space of square-integrable complex psi-functions, i.e., the Hilbert 
space -L^(^). 

When the system is in the state \ip{t)), the average value (/) of any physical observable 
/ is equal to 

(/) = m)\fm)). 

where / is the Hermitian operator corresponding to /. 

A quantum system is coherent if it is in a linear superposition of its basis states. If a 
measurement is performed on the system and this means that the system must somehow 
interact with its environment, the superposition is destroyed and the system is observed to 
be in only one basis state, as required classically. This process is called reduction or collapse 
of the wavefunction or simply decoherence and is governed by the form of the wavefunction 

Entanglement, on the other hand, is a purely quantum phenomenon and has no classical 
analogue. It accounts for the ability of quantum systems to exhibit correlations in counter- 
intuitive 'action-at-a-distance' ways. Entanglement is what makes all the difference in the 
operation of quantum computers versus classical ones. Entanglement gives 'special powers' 
to quantum computers because it gives quantum states the potential to exhibit and maintain 
correlations that cannot be accounted for classically. Correlations between bits are what 
make information encoding possible in classical computers. For instance, we can require two 
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bits to have the same value thus encoding a relationship. If we are to subsequently change the 
encoded information, we must change the correlated bits in tandem by explicitly accessing 
each bit. Since quantum bits exist as superpositions, correlations between them also exist in 
superposition. When the superposition is destroyed (e.g., one qubit is measured), the correct 
correlations are instantaneously 'communicated' between the qubits and this communication 
allows many qubits to be accessed at once, preserving their correlations, something that is 
absolutely impossible classically. 

More precisely, the first quantization is a linear representation of all classical dynamical 
variables (like coordinate, momentum, energy, or angular momentum) by linear Hermitian 
operators acting on the associated Hilbert state-space 7i, which has the following properties 

m- 

1. Linearity: af + (3g ^ a f + (3 g, (for all constants a, (3 £ C); 

2. A 'dynamical' variable, equal to unity everywhere in the phase-space, corresponds to 
unit operator: 1^1; and 

3. Classical Poisson brackets 

df dg df dg 



if ^9} Qqi dpi dpidq^ 
quantize to the corresponding commutators 

{/, 9} ^ g] , [f, 9] = fg- gf- 

Like Poisson bracket, commutator is bilinear and skew-symmetric operation, satisfying 
Jacobi identity. For Hermitian operators f,g their commutator [f,g] is anti-Hermitian; for 
this reason i is required in {/,<?} — > —ih[f,g]. 

Property (2) is introduced for the following reason. In Hamiltonian mechanics each dy- 
namical variable / generates some transformations in the phase-space via Poisson brackets. 
In quantum mechanics it generates transformations in the state-space by direct application 
to a state, i.e., 

u = {u,f}, dt\^) = ^m- (2) 

Exponent of anti-Hermitian operator is unitary. Due to this fact, transformations, gen- 
erated by Hermitian operators 

U = exp— , 

are unitary. They are motions - scalar product preserving transformations in the Hilbert 
state-space Ti. For this property i is needed in ([2]). 
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Due to property (2), the transformations, generated by classical variables and quantum 
operators, have the same algebra. 

For example, the quantization of energy E gives: 

E ^ E = ihdt. 

The relations between operators must be similar to the relations between the relevant physical 
quantities observed in classical mechanics. 

For example, the quantization of the classical equation E = H, where 

H = H{p„q')=T + U 

denotes the Hamilton's function of the total system energy (the sum of the kinetic energy 
T and potential energy U), gives the Schrodinger equation of motion of the state ket-vector 
\ip{t)) in the Hilbert state-space Ti. 

ihdtim) =H\m)- 

In the simplest case of a single particle in the potential field C/, the operator of the total 
system energy - Hamiltonian is given by: 

Zm 

where m denotes the mass of the particle and V is the classical gradient operator. So the 
first term on the r.h.s denotes the kinetic energy of the system, and therefore the momentum 
operator must be given by: 

p = —\hV. 

Now, for each pair of states \ip), their scalar product is introduced, which is [IB] : 

1. Linear (for right multiplier): 

((^laiV-i + "2V'2) = ai(¥'IV'i) + a2(<^|V'2); 

2. In transposition transforms to complex conjugated: 



{Lp\ip) = {il;\ip)- 
this implies that it is 'anti-linear' for left multiplier: 
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3. Additionally it is often required, that the scalar product should be positively defined: 

forahlV'), {ip\i^)>0 and (ipi^) = iff |V) = 0. 

Complex conjugation of classical variables is represented as Hermitian conjugation of 
operators!! 

If the two Hermitian operators / and g commute, i.e., [f,g] = (see Heisenberg picture 
below), than the corresponding quantities can simultaneously have definite values. If the two 
operators do not commute, i.e., [f,g] ^ 0, the quantities corresponding to these operators 
cannot have definite values simultaneously, i.e., the general Heisenberg uncertainty relation 
is valid: 

{Af)'.{Agf>^^[f,gf, 

where A denotes the deviation of an individual measurement from the mean value of the dis- 
tribution. The well-known particular cases are ordinary uncertainty relations for coordinate- 
momentum {q — p), and energy-time {E — t): 

Aq-Apq>^, and AE-At>^. 

For example, the rules of commutation, analogous to the classical ones written by the 
Poisson's brackets, are postulated for canonically-conjugate coordinate and momentum op- 
erators: 

^Two operators /, are called Hermitian conjugated (or adjoint), if 

Mfi>)^{f^vW (forall<^,V). 

This scalar product is also denoted by {if\f\tp) and called a matrix element of an operator. 

- operator is Hermitian (self-adjoint) if = / and anti-Hermitian if /"'" = — /; 

- operator is unitary, if = U~^; such operators preserve the scalar product: 

Real classical variables should be represented by Hermitian operators; complex conjugated classical variables 
(a, a) correspond to Hermitian conjugated operators (a,a+). 

Multiplication of a state by complex numbers does not change the state physically. 
Any Hermitian operator in Hilbert space has only real eigenvalues; 

M) = MA), (for all /. GK). 

Eigenvectors \ form complete orthonormal basis (eigenvectors with different eigenvalues are automatically 
orthogonal; in the case of multiple eigenvalues one can form orthogonal combinations; then they can be 
normalized). 
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where 5*- is the Kronecker's symbol. By applying the commutation rules to the system 
Hamiltonian H = H{pi,(f'), the quantum Hamilton's equations are obtained: 

Pi = -dqrH and = dp.H. 

A quantum state can be observed either in the coordinate q— representation, or in the 
momentum p— representation. In the g— representation, operators of coordinate and momen- 
tum have respective forms: q = q, and pq = —ih-^, while in the p-representation, they have 

respective forms: q = i^g^, and Pq = Pq. The forms of the state vector \ip{t)) in these two 
representations are mathematically related by a Fourier-transform pair. 

1.1.2 Three quantum pictures 

In the representation, there are three main pictures (see e.g., |16)): 

1. Schrddinger picture, 

2. Heisenberg picture, and 

3. Dirac interaction picture. 

These three pictures mutually differ in the time-dependence, i.e., time-evolution of the 
state-vector wavefunction \ip{t)) and the Hilbert coordinate basis (g*) together with the 
system operators. 

1. In the Schrddinger (S) picture, under the action of the evolution operator S{t) the 
state-vector \ip{t)) rotates: 

\m) =s{t)\m), 

and the coordinate basis (g*) is fixed, so the operators are constant in time: 

F{t) = F(0) = F, 

and the system evolution is determined by the Schrodinger wave equation: 

ihdt\iJ^'^{t)) = H^\i;^'^{t)). 

If the Hamiltonian does not explicitly depend on time, H{t) = H, which is the case with 
the absence of variables of macroscopic fields, the state vector \ip{t)) can be presented in the 
form: 

IV'(t)) =exp(-iff) \tl;), 
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satisfying the time-independent Schrodinger equation 



H\i^) =E\il,), 

which gives the eigenvalues Em and eigenfunctions of the Hamiltonian H. 

2. In the Heisenberg (H) picture, under the action of the evolution operator S(t), the 
coordinate basis (g*) rotates, so the operators of physical variables evolve in time by the 
similarity transformation: 

F{t) = S-\t)F{0)S{t), 
while the state vector \tp{t)) is constant in time: 

\m) = \m) = m, 

and the system evolution is determined by the Heisenberg equation of motion: 

ihdtF''it) = [F"it),H''it)], 

where F{t) denotes arbitrary Hermitian operator of the system, while the commutator, i.e., 
Poisson quantum bracket, is given by: 

[F{t), H{t)] = F{t) H{t) - H{t) F{t) = iK. 

In both Schrodinger and Heisenberg picture the evolution operator S{t) itself is determined 
by the Schrodinger-like equation: 

ihdtSit) = HS{t), 

with the initial condition ^(O) = /. It determines the Lie group of transformations of the 
Hilbert space L^(^) in itself, the Hamiltonian of the system being the generator of the group. 

3. In the Dirac interaction (I) picture both the state vector \ip{t)) and coordinate basis (g*) 
rotate; therefore the system evolution is determined by both the Schrodinger wave equation 
and the Heisenberg equation of motion: 

ihdt\2p'{t)) = \^\t)), and \hdtF\t) = [F\t),HO{t)]. 

Here, H = + , where corresponds to the Hamiltonian of the free fields and 
corresponds to the Hamiltonian of the interaction. 

In particular, the stationary (time-independent) Schrodinger equation. 

Hip = Eip, 

can be obtained from the condition for the minimum of the quantum action: 

SS[tlj] = 0. 
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The quantum action is usually defined by the integral: 



S[i;]= {m\H\m) = jrH^ dV, 
with the additional normalization condition for the unit-probability of the psi-function: 

(mm)) = I r^dv = i. 



When the functions if) and ip* are considered to be formally independent and only one of 
them, say is varied, we can write the condition for an extreme of the action: 

SS[^] = [ 5il}*Hi) dV-E [ 6ij*ij dV = [ 5^*{H^ - Ei;) dV = 0, 



where E' is a Lagrangian multiplier. Owing to the arbitrariness of (5V'*) the Schrodinger 
equation ~ -^V' = must hold. 



1.1.3 Dirac's probability amplitude and time dependent perturbation 

Most quantum-mechanical problems cannot be solved exactly. For such problems we can use 
Dirac's perturbation method, which consists in splitting up the time-dependent Hamiltonian 
H = H{t) into two parts: 

H{t) = HQ + eHi{t), 

in which Hq = E must be simple, non-autonomous, energy function that can be dealt with 
exactly, while eHi{t) = V{t) is small time-dependent perturbation, which can be expanded 
as a power series in a small numerical factor e. The first part, Hq, may then be consid- 
ered as the Hamiltonian of a simplified, or unperturbed system that can be exactly solved, 
while the addition of the second part eHi{t) will require small corrections, of the nature 
of a power-series expanded perturbation in the solution for the unperturbed system. Pro- 
vided the perturbation series in e converges, the perturbation method will give the answer 
to our problem with any desired accuracy; even when the series does not converge, the first 
approximation obtained by means of it is usually fairly accurate [15] . 

Therefore, we do not consider any modification to be made in the states of the unperturbed 
system E = Hq, but we suppose that the perturbed system H{t), instead of remaining 
permanently in one of its states, is continually changing from one state to another (or, making 
transmissions), under the influence of the perturbation V{t) = eHi{t). 

We will work in the Heisenberg's representation for the unperturbed system E, assuming 
that we have a general set of linear Hermitian operators a's to label the representatives. Let 
us suppose that at initial time to the system is in a state for which the a's certainly have 
the values a', so that the basic ket \a') would correspond to this state. This state would be 
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stationary if there were no perturbation, i.e., if H{t) = E. The perturbation V{t) cause the 
E to change. At time t the ket corresponding to the state \a') in the Schrddinger's picture 
will be T I a') , according to equation 

\At) = T \ Ato) , as well as 

i^ = H{t)T and i^l^ = H{t) \At) , (3) 

where T is a linear Hermitian operator independent of the ket \At) and depending only on 
time (to and t). The probability of the a's having the values a" is given by the absolute 
square of the probability amplitude (or, transition amplitude) {a"\T\a') (for the system's 
transition from the state \a') to the state \a")), 

P{a',a") = \{a"\T\a')\\ (4) 

For a' 7^ a", P{a',a") is the probability of a transition taking place from state a' to state 
a" during the time interval [toi j P{oi\ is the probability of no transition taking place at 
all, while the sum of P{a' , a") for all a" is unity. 

Let us now suppose more generally that initially the system is in one of the various states 
a' with the probability Pa' for each. To deal effectively with this problem, we introduce 
the von Neumann's quantum density function p, a quantum-mechanical analogue to the 
Gibbs statistical density function p = p[t) of a Gibbs ensemble with the classical Hamiltonian 
H(q,p), which evolves within the ensemble's n— dimensional phase-space V = {{q^,Pi) \ i = 
l,...,n} according to the Poisson equation 

dtp = -[p,H{q,p)], 
with the normalizing condition : J J pdq^dpi = 1. 

The von Neumann's quantum density function p corresponding to the initial probability 
distribution P^' is given by 

Po = '^\(^') Pa' {a\- 

a' 

At time t, each ket \a') will have changed to T\a') and each bra {a'\ will change to {a'\ T 
(where T is complex-conjugate to T), so p^ will have changed to 

p{t) = Y,T\a')Pa'{a'\f. 

a' 

The probability amplitude of a's then having the values of a" will be (using ^) 
{a"\ p{t) \a') = Y, T \a') P^' {a'\ f \a") = Pa'P{a', a"). 
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This result expresses that the probabihty of the system being in the state a" at time t equals 
the sum of the probabilities of the system being initially in any state a' ^ a" , and making 
a transition from state a' to the final state a" . Thus, the various transition probabilities act 
independently of one another, according to the ordinary laws of probability. 

The whole problem of calculating transitions thus reduces to the determination of the 
probability amplitudes {a"\T\a') [15j. These can be worked out from ([3]), or 

if = [Ho + eHi{t)]T = {E + V)T (where t = dT/dt). (5) 

This calculation can be simplified if instead of T and V operators, we are working with 

T* = exp[i£;(t - tQ)]T and V* = exp[i^(t - to)]V eiq)[-iE{t - to)], (6) 

where V* is the result of applying a unitary transformation to V. Using ([6]) we obtain 

if* =e^^[iE{t-tQ)]VT = V*T\ (7) 

which is more convenient then ^ as it makes the time evolution of T* depend only on the 
(unitary transformed) perturbation V* and not on the unperturbed state E. From ([U]) we get 
the probability amplitude 

{a"\ T* \a') = exp[i£;(t - to)] (a"| T\a') , so that 
P{a\a") = |(a"|r*|a'>|2, 

which shows that T and T* are equally good for determining transition probabilities. 

So far, our work in this subsection has been exact. Now we assume the perturbation 
V{t) = eHi{t) is a small quantity of the first order in e and express T* in the form 

T* = 1 + Ti* +r2* + (8) 

where = (e), = T^e^), etc. Substituting ^ into ^ we get the expansion 

if I = V*, 
if* = V*T^, 
if* = V*T^, 

From the first of these equations we obtain 

= -i [ V*{t')dt', 

J to 



14 



from the second we get 
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V*{t')dt' [ V*{t")dt 



II 



and so on. 



Now, the perturbation form of the transition probabihty P(a',a") = | (q"| T* \ol) ^ is, 
if we retain only the first-order term (which is sufficiently accurate for many practical 
problems), given by 



If we retain the first two terms, Tj* and T2 , the transition probability -P(a', a") is given by 



where a'" is the so-called intermediate state (between a' and a"). This shows the perturbative 
calculation of the transition probability P{a' , a"): we fist calculate the perturbative expansion 
of the transition amplitude {a"\ T* \ a') , and then take its absolute square to obtain the overall 
transition probability. For more technical details, see p!5] . 

Both Dirac's concepts introduced in this subsection, namely transition amplitude and 
time— dependent perturbation, will prove essential later in the development of the Feynman 
path integral, as well as the Feynman diagrams approach to quantum field theory (QFT). 

1.1.4 State space for n non-relativistic quantum particles 

Classical state-space for the system of n particles is its 6iVD phase-space V, including all 
position and momentum vectors, rj = {x,y,z)i and Pi = {px,Py,Pz)i respectively (for i = 
1, ...,n). The quantization is performed as a linear representation of the real Lie algebra Cp 
of the phase-space V, defined by the Poisson bracket {A, B} of classical variables A,B ~ into 
the corresponding real Lie algebra Ch of the Hilbert space 7i, defined by the commutator 
[A, B] of skew-Hermitian operators A, B [16j . 

We start with the Hilbert space TCx for a single ID quantum particle, which is composed 
of all vectors of the form 






ip (x) |x) dx, 



where (x) = {x\ip) are square integrable Fourier coefficients. 
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The position and momentum Hermitian operators, x and p, respectively, act on the vectors 
IV'a;) S T^x in the following way: 



/ + 00 /• + 00 

X ip {x) \x) dx , I \x il) {x)\ dx < +00 , 
-00 J —00 

/ + 00 f + OD 

—ihdx'4' (x) \x) dx, / \—ihdx'(p {x)\^ dx < +00. 
-00 J —00 

The orhit Hilbert space TC^ for a single 3D quantum particle with the full set of compatible 
observable f ={x,y, z),p = {px,Py,Pz), is defined as 

where f has the common generalized eigenvectors of the form 

|f) = \x)x\y)x\z) . 
Hi is composed of all vectors l'^^) of the form 

/• z'+oo /»4-oo Z'+OO 

\ipr) = / V'(r) \r) dr= / / il){x,y,z) \x)x\y)x\z) dxdydz, 

where ip (r) = (rl'^^) are square integrable Fourier coefficients. 



/+00 
|^(r)| ^dr < +00. 
-00 



The position and momentum operators, f and p, respectively, act on the vectors £ ^1 
in the following way: 



= / f ■?/) (r) |r) dr, / |r (r)| ^ dr < +00, 
Jn° JHl 

IV'^) = / -i/iSf V (r) |r) dr, / {r)\^ dr < +00 . 



Now, if we have a system of n 3D particles, let HI denote the orbit Hilbert space of the 
ith particle. Then the composite orbit state-space 7^° of the whole system is defined as a 
direct product 

7^° is composed of all vectors 

\i>r)= / 'i/'(ri,r2, ...,r„) |ri)x|r2)x...x|r„) dridr2...dr„ 
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where (ri,r2, ...,r„) = (ri,r2, ...,r„|'0") are square integrable Fourier coefficients 

/ |V'(ri,r2,...,r„)|^ dridr2...drn < +00. 

The position and momentum operators and Pi act on the vectors \ipl^) G H!^ in the following 
way: 

^Mr) = / {fi}'0(ri,r2,...,r„) |ri)x|r2)x...x|r„) dridr2...dr„, 
pMv) = / {-i/i5fi}V'(ri,r2,...,r„) |ri)x|r2)x...x|r„) dri(ir2...dr„, 
with the square integrable Fourier coefficients 

/ |{fi}V'(ri,r2,...,r„)|^ dridr2...dr„ < +00, 

/ \{-ihdri}ip{ri,r2,...,rn)f dridr2...drn < +00, 

respectively. In general, any set of vector Hermitian operators {Aj} corresponding to all the 
particles, act on the vectors £ in the following way: 

AilV'r) = / {Aij-^ (ri,r2, ...,r„) |ri)x|r2)x...x|r„) dridr2...dr„, 
with the square integrable Fourier coefficients 



/ |{Ai} V'(ri,r2,...,i 



2 

dridr2...drn < +00. 



1.2 Transition to quantum fields 

1.2.1 Amplitude, relativistic invariance and causality 

We will see later that in QFT, the fundamental quantity is not any more Schrodinger's 
wavefunction but the rather the closely related, yet different, Dirac-Feynman's amplitude. 
To introduce the amplitude concept within the non-relativistic quantum mechanics, suppose 
that X = (a:;, y, z) and consider the amplitude for a free particle to propagate in time t from 
xq to X = x(t), which is given by 

[/(t) = <x|e-'^*|xo). 
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As the kinetic energy of a free particle is E = p^/2m, we have 



U{t) = (x|e-^(P'/2-)*|xo) = / d3p/x|e-'(p'/2'")*|p)(p|xo) 



= ^ -i(p2/2m)t ip-(x-xo) ^ im(x-xo)2/2t 

(27r)2 y ^ \27ritJ 

Later we will deal with the amplitude in the relativistic framework. 

As we have seen in the previous section, in non-relativistic quantum mechanics observables 
are represented by self-adjoint operators that in the Heisenberg picture depend on time, so 
that for any quantum observable O, the Heisenberg equation of motion holds (in normal 
units): 

idtO = [0,H]. (9) 

Therefore measurements are localized in time but are global in space. The situation is 
radically different in the relativistic case. Because no signal can propagate faster than the 
speed of light, measurements have to be localized both in time and space. Causality demands 
then that two measurements carried out in causally-disconnected regions of space-time cannot 
interfere with each other. In mathematical terms this means that if Or-^ and O/jj are the 
observables associated with two measurements localized in two causally-disconnected regions 
Ri, R2, they satisfy the commutator relation 

[Ob,,,Or2]=0, if (xi - X2)^ < 0, for all xi G X2 G i?2. (10) 

Hence, in a relativistic theory, the basic operators in the Heisenberg picture must depend 
on the space-time position x^. Unlike the case in non-relativistic quantum mechanics, here 
the position x is not an observable, but just a label, similarly to the case of time in ordinary 
quantum mechanics. Causality is then imposed microscopically by requiring 

[Oix),Oiy)]=0, if(x-y)2<0. (11) 

A smeared operator Or over a space-time region R can then be defined as 

Or = j d^xO{x)fR{x), 

where fR{x) is the characteristic function associated with R, 



fdx) 



1, for X £ R, 
0, for X ^ R. 



Equation (jlOp follows now from the micro-causality condition (jlip . 

Therefore, relativistic invariance forces the introduction of quantum fields. It is only when 
we insist in keeping a single-particle interpretation that we crash against causality violations. 
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To illustrate the point, let us consider a single particle wave function ip{t,x) that initially is 
localized in the position x = 

^'(O, x) = 5{x). 



Evolving this wave function using the Hamiltonian H = \/— + m? we find that the 
wave function can be written as 

Integrating over the angular variables, the wave function can be recast in the form 

The resulting integral can be evaluated using the complex integration contour C. The result 
is that, for any f > 0, one finds that 'ip{t,x) ^ for any x. If we insist in interpreting the 
wave function x) as the probability density of finding the particle at the location x in 
the time t we find that the probability leaks out of the light cone, thus violating causality. 

In the Heisenberg picture, the amplitude for a particle to propagate from point y to point 
X in the field (p is defined as 

D{x-y) = {0\<P{x)Hym = J ^^e-H-^/). 

In this picture we can make time-dependent the field operator (j) = <t>{x) and its cannonically- 
conjugate momentum operator vr = '/r(x), as 

(t){x) = ,/.(x, t) = ei^V(x) 'Kix) = 7r(x, t) = e'-^*7r(x) e'^^*. 

Using ([9]) we can compute the time dependence of 4> and vr as 

i0(x, t) = i7r(x, t), i7r(x, t) = -i(-V^ + m^)(/)(x, t). 

Combining the two results we get the Klein-Gordon equation 

ij)= (v2-m2)(/). 



1.2.2 Gauge theories 

Recall that a gauge theory is a theory that admits a symmetry with a local parameter. For 
example, in every quantum theory the global phase of the wave ^—function is arbitrary and 
does not represent something physical. Consequently, the theory is invariant under a global 
change of phases (adding a constant to the phase of all wave functions, everywhere); this is 
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a global symmetry. In quantum electrodynamics, the theory is also invariant under a local 
change of phase, that is, one may shift the phase of all wave functions so that the shift may 
be different at every point in space-time. This is a local symmetry. However, in order for 
a well-defined derivative operator to exist, one must introduce a new field, the gauge field, 
which also transforms in order for the local change of variables (the phase in our example) not 
to affect the derivative. In quantum electrodynamics this gauge field is the electromagnetic 
potential 1-form A (see Appendix), in components within the nD coframe {dx^} on a smooth 
manifold M (dual to the frame, i.e., basis of tangent vectors {d^ = d/dx^}, given by 

A = A^dx^, such that ^new = ^oid + df, (/ is any scalar function) 

- leaves the electromagnetic field 2-form F = dA unchanged. This change df of local gauge 
of variable A is termed gauge transformation. In quantum field theory the excitations of 
fields represent particles. The particle associated with excitations of the gauge field is the 
gauge boson. All the fundamental interactions in nature are described by gauge theories. 
In particular, in quantum electrodynamics, whose gauge transformation is a local change 
of phase, the gauge group is the circle group U{1) (consisting of all complex numbers with 
absolute value 1), and the gauge boson is the photon (see e.g., [T7]). 

The gauge field of classical electrodynamics, given in local coframe {dx^} on M as an 
electromagnetic potential 1-form 

A = A^dx^ = A^dx^ + df, (/ = arbitrary scalar field), 

is globally a connection on a ;7(l)-bundle of M0 The corresponding electromagnetic field, 

''Recall that in the 19th Century, Maxwell unified Faraday's electric and magnetic fields. Maxwell's theory 
led to Einstein's special relativity where this unification becomes a spin-off of the unification of space end time 
in the form of the Faraday tensor [18j 

F = E Adt + B, 

where F is electromagnetic 2— form on space-time, E is electric 1— form on space, and B is magnetic 2— form 
on space. Gauge theory considers F as secondary object to a connection-potential 1— form A. This makes 
half of the Maxwell equations into tautologies P2], i.e., 

F = dA dF ^ : Bianchi identity, 

but does not imply the second half of Maxwell's equations, 

5F — — IttJ : dual Bianchi identity. 
To understand the deeper meaning of the connection-potential 1— form A, we can integrate it along a path 7 

in space-time, x *- y. Classically, the integral A represents an action for a charged point particle to 

move along the path 7. Quantum-mechanically, exp ^i AlJ represents a phase (within the unitary Lie group 

f/(l)) by which the particle's wave-function changes as it moves along the path 7, so ^4 is a [7(1)— connection. 

In other words, Maxwell's equations can be formulated using complex line bundles, or principal bundles with 
fibre (7(1). The connection V on the line bundle has a curvature F = V'^ which is a 2-form that automatically 
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locally the 2-form on M, 



F 
F 



dA, 
1 



in components given by 



Fr,^ dx^" A dx" , with F, 



is globahy the curvature of the connection under the gauge-covariant derivative, 

D^ = d^-ieA^, (12) 

where e is the charge couphng constant]! In particular, in 4D space-time electrodynamics, 
the 1-form electric current density J has the components = = (p, jx, jy, jz) (where 
p is the charge density) , the 2-form Faraday F is given in components of electric field E and 
magnetic field B by 



^ IIV 



/ 





Ex 


Ey 


E, 


\ 











By 






-Ey 


Bz 





—Bx 




V 




-By 


Bx 





/ 



with 



-F 



satisfies dF = and can be interpreted as a field-strength. If tfie line bundle is trivial with fiat reference 
connection d, we can write V = d + A and F = dA with A the 1-form composed of the electric potential and 
the magnetic vector potential. 

^The only thing that matters here is the difference a between two paths and 72 in the action J A [19) . 
which is a 2-morphism (see e.g., [201 121) ) 




®If a gauge transformation is given by 
and for the gauge potential 

then the gauge-covariant derivative, 
transforms as 
and D^tp transforms as 



Au 



A^ + -{d^A), 

e 



= df_,- ieA^ 



Df_,^ dfj.~ ieA^ - i(9pA). 
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while its dual 2-form Maxwell -kF has the following components 



( —Bx —By -Bz \ 

Bj. —Ez Ey 

By Ez —Ex 

\ Bz —Ey Bx / 



with Fy^ = - -k F^y, 



so that classical electrodynamics is governed by the Maxwell equations^ which in modern 
exterior formulation read 

dF = 0, 5F = —4:71 J, or in components, 

where * is the Hodge star operator and 5 is the Hodge codiferential (see section 13.21 be- 
low), comma denotes the partial derivative and the 1-form of electric current J = J^dx^ is 
conserved, by the electrical continuity equation, 

5 J = 0, or in components, J^,'^ = 0. 

The first, sourceless Maxwell equation, dF = 0, gives vector magnetostatics and magne- 
todynamics. 

Magnetic Gauss' law : divB = 0, 

Faraday's law : S^B + curlE = 0. 

The second Maxwell equation with source, 6F = J, gives vector electrostatics and electrody- 
namics, 

Electric Gauss' law : divE = 47rp, 

Ampere's law : 5jE — curlB = — 47rj. 

The standard Lagrangian for the free electromagnetic field, F = dA, is given by |20 ^ l21 t [T6] 

£(A) = i(FA*F), 
with the corresponding action functional 

S{A) = ^ I FA^F 

Maxwell's equations are generally applied to macroscopic averages of the fields, which 
vary wildly on a microscopic scale in the vicinity of individual atoms, where they undergo 
quantum effects as well (see below). 



22 



1.2.3 Free and interacting field theories 

A generic gauge-covariant derivative with Lorentz index /x is denoted by Z)^. For a Maxwell 
field, Dfj, is given by p^ . 

Dirac slash notation : ^ =^ 1^9^, p =^ l^D^, 

Dirac algebra : {7^^, 7^} = 7'^7^ + 7'^7'^ = 2^'^^ x l^xn- 

Standard free theories are Klein-Gordon and Dirac fields: 

Klein-Gordon equation : (9^ + m'^)ilj = 0, 
Dirac equation : (i^ — m)i/^ = 0. 

Two main examples of interacting theories are (/)^— theory and QED: 

1. (/>^— theory: 

Lagrangian : C = ^{df,(j)f - ^m^(j)^ - ^cp^, 
Equation of motion : (9^ + m'^)(j) = 

2. QED: 

Lagrangian : £ = ^Maxwell + ^Coirac + J^int 

= -^iF^,f + i;iiip-m)^p. 

Gauge invariance : ipix) — > e'^^^VC^^) =^ ~^ dna{x). 

e 

Equation of motion : (i-^— m)'ip = 0. 

1.2.4 Dirac QED 

The Dirac equation for a particle with mass m (in natural units) reads (see, e.g., [16j ) 

(17/^9^ - m)V^ = 0, (/u = 0,l,2,3) (13) 
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where tlj{x) is a 4-component spinoitj wave-function, the so-called Dirac spinor, while 7^ are 
4x4 Dirac j— matrices, 
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They obey the anticommutation relations 

where g^^ is the metric tensor. 

Dirac's 7— matrices are conventionally derived as 

where are Pauli a—matrice^ (a set of 2 x 2 complex Hermitian and unitary matrices), 
defined as 

obeying both the commutation and anticommutation relations 

[cri,aj] = TiEijkCTk, {ai,aj} = 25ij ■ I, 

where Eijk is the Levi-Civita symbol, 6ij is the Kronecker delta, and / is the identity matrix. 

Now, the Lorentz-invariant form of the Dirac equation (jl3p for an electron with a charge 
e and mass me, moving with a 4-momentum 1-form p = p^dx^ in a classical electromagnetic 
field defined by 1-form A = A^dx'^, reads (see, e.g., [T6]): 

{i7^ [p^ - eA^\ - me} = 0, (14) 
'^The most convenient definitions for the 2-spinors, like the Dirac spinor, are: 



"1 




,0^ = 


1 


and — 


O' 
1 




l' 




*In quantum mechanics, each Pauli matrix represents an observable describing the spin of a spin i particle 
in the three spatial directions. Also, ia^ are the generators of rotation acting on non-relativistic particles with 
spin I . The state of the particles are represented as two-component spinors. 

In quantum information, single-qubit quantum gates are 2x2 unitary matrices. The Pauli matrices are 
some of the most important single-qubit operations. 
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and is called the covariant Dirac equation. 

The formal QED Lagrangian (density) includes three terms, 

C{x) = Cemix) + Cint{x) + Cc-p{x), (15) 

related respectively to the free electromagnetic field 2-form F = F^i,dx^ Adx'^ , the electron- 
positron field (in the presence of the external vector potential 1-form A^^^), and the interac- 
tion field (dependent on the charge-current 1-form J = J^dx^). The free electromagnetic 
field Lagrangian in (jl5p has the standard electrodynamic form 

r (r) — --F^^F 

where the electromagnetic fields are expressible in terms of components of the potential 1- 
form A = A^dx^ by 

F^, = d^A'j"-d,Al°\ with 4°* = ^-* + ^^. 

The electron-positron field Lagrangian is given by Dirac's equation (|14p as 

£e-p(x) = i;{x) {iT'^ - e^^*] - m,} ^{x), 

where 'ip{x) is the Dirac adjoint spinor wave function. 
The interaction field Lagrangian 

accounts for the interaction between the uncoupled electrons and the radiation field. 
The field equations deduced from (fT5|) read 

{iT'^b^-eyl^^*] -me}V(x) = l^Hx)A^, 

d^F^, = J,. (16) 

The formal QED requires the solution of the system p6|) when A^{x), ij){x) and 'tp{x) are 
quantized fields. 



1.2.5 Abelian Higgs Model 

The Abeliarll Higgs model is an example of gauge theory used in particle and condensed 
matter physics. Besides the electromagnetic field it contains a self-interacting scalar field, 

®An Abelian (or, commutative) group (even better, Lie group, see Appendix), is such a group G that 
satisfies the condition: a ■ b — b ■ a for ail a,b £ G. In other words, its commutator, [a, b] :— a^^b^^ab equals 
the identity element. 
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the so-cahed Higgs field, minimally coupled to electromagnetism. From the conceptual point 
of view, it is advantageous to consider this field theory in (2 + 1)D space-time and to extend 
it subsequently to (3 + 1)D for applications. The Abelian Higgs Lagrangian reads [39] 

which contains the complex (charged), self-interacting scalar field (p. The Higgs potential is 
the Mexican hat function of the real and imaginary part of the Higgs field. 

By construction, this Higgs potential is minimal along a circle \(p\ = a in the complex (p 
plane. The constant A controls the strength of the self-interaction of the Higgs field and, for 
stability reasons, is assumed to be positive, A > . The Higgs field is minimally coupled to 
the radiation field A^, i.e., the partial derivative 9^ is replaced by the covariant derivative, 
DfM = dfj_ + ieA^ . Gauge fields and field strengths are related by 

F^,, = d^A, - d,A^ = - [D^, D,] . 

le 

The inhomogeneous Maxwell equations are obtained from the least action principle, 

6S = 0, with S = [ Cd^x = 0, 



by variation of the action S with respect to the gauge fields A^ (and their derivatives df^A^). 
With 

-F^^ TT- = -3 ' we get 



5d^A, ' 5A, 



dpF'^" = f, 3u = H<^'d,<P-(t>d,(t>^)-2e'<P*<pA,. 

We remark here that the homogeneous Maxwell equations are not dynamical equations 
of motion. They are integrability conditions and guarantee that the field strength can be 
expressed in terms of the gauge fields. The homogeneous equations follow from the Jacobi 
identity of the covariant derivative 

[D^, [D,, D^W + [D„, [D^, D,]] + [Z?„ [D„, D^]] = 0. 

Multiplication with the totally antisymmetric 4-tensor e^^'^^ , yields the homogeneous equa- 
tions for the dual field strength F^^ 



0, F^^^ = ^^e^'^P^Fp^. 
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The transition: F — > F corresponds to the following duality relation of electric and 
magnetic fields, E — > B, B — > — E. 

Variation with respect to the charged matter field yields the equation of motion: 

«„«.^ + |^ = o. 

Gauge theories contain redundant variables. This redundancy manifests itself in the 
presence of local symmetry transformations, or gauge transformations, U{x) = e^^°'^^\ which 
rotate the phase of the matter field and shift the value of the gauge field in a space-time 
dependent manner 

^^^[U] = u{x)m , aP =A^ + U{x) i U\x) . (17) 

The covariant derivative -D^ has been defined such that Dfj^cp transforms covariantly, i.e., like 
the matter field (f) itself. 

D^(j){x) -^U{x)D^ct){x). 

This transformation property together with the invariance of F^i, guarantees invariance of C 
and of the equations of motion. A gauge field which is gauge equivalent to A^ = is called 
a pure gauge. According to (fTTj) a pure gauge satisfies 

Affix) = U{x)^ U^{x) = -d^ a(x), 

and the corresponding field strength vanishes. 

Note that the non-Abelian Higgs model has the action: 

S{<I),A) = ^ I TT{F>^''F,,) + \Dcl>\'' + vm), 

where now the non-Abelian field A is contained both in the covariant derivative D and in the 
components F^^'^ and Fjj^i^ (see Yang-Mills theory below). 

2 Feynman Path Integral 

2.1 The action— amplitude formalism 

The 'driving engine' of quantum field theory is the Feynman path integral. Very briefiy, there 
are three basic forms of the path integral (see, e.g., |22l I16j ): 
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1. Sum-over-histories, developed in Feynman's version of quantum mechanics (QM0 

2. Sum-over-fields, started in Feynman's version of quantum electrodynamics (QED) 

isaj 

and later improved by Fadeev-Popov [25] : 

3. Sum-over-geometries/topologies in quantum gravity (QG), initiated by S. Hawking 
and properly developed in the form of causal dynamical triangulations (see for a 'softer' 
review, see [27]). 

In all three versions, Feynman's action- amplitude formalism includes two components: 



^"Feynman's amplitude is a space-time version of the Schrodinger's wavefunction ip, which describes how 
the (non-relativistic) quantum state of a physical system changes in space and time, i.e., 

{Outt^;^|Int;^J = i/'(x,t), (for X G [In, Out], t G [ti„i,t/i„]). 

In particular, quantum wavefunction tp is a complex-valued function of real space variables x — 
{xi, X2, ^n) G R", which means that its domain is in R" and its range is in the complex plane, formally 
i/;(x) : R" —> C. For example, the one-dimensional stationary plane wave with wave number k is defined as 

tpix) = e"""", (for X G R), 

where the real number k describes the wavelength, \ = 2n /k. In n dimensions, this becomes 

V(:E)=e'P■^ 

where the momentum vector p = k is the vector of the wave numbers k in natural units (in which h — m — 1). 

More generally, quantum wavefunction is also time dependent, ip = i/)(x,t). The time-dependent plane wave 
is defined by 

i/;(x,t) =e'P'"-'f'*/^ (18) 
In general, ^'(''^i is governed by the Schrodinger equation |28l 116) (in natural units /i = m = 0) 

i^V(x,t) = -iAV'(x,t), (19) 

where A is the n— dimensional Laplacian. The solution of (|19|l is given by the integral of the time-dependent 
plane wave ifTSjl. 

^(x,t) = ^£^e--^^*/^^o(p)d>, 
which means that t) is the inverse Fourier transform of the function 

V'(p,t)=e-'^'*''^^o(p), 

where </'o(p) calculated for each initial wavefunction. For example, if initial wavefunction is Gaussian, 

2;2 - 1 

/(x) = exp(— a — ), with the Fourier transform /(p) = — = exp( ). 

2 ■ Ja 2a 



then ^o(p) = exp(--f^) 
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1. A real-valued, classical, Hamilton's action functional, 

:= / L[^]dt, (20) 
with the Lagrangian energy function defined over the Lagrangian density C, 

while $ is a common symbol denoting all three things to be summed upon (histories, fields 
and geometries). The action functional S[^] obeys the Hamilton's least action principle, 
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5S'[$] = 0, and gives, using standard variational methods^ the Euler-Lagrangian equations, 
which define the shortest path, the extreme field, and the geometry of minimal curvature (and 
without holes). 



In Lagrangian field theory, the fundamental quantity is the action 



S[$] = f Ldt^ I d"xJZ{$,d^$) 
Jti„ Jm.^ 



so that the least action principle, ^Sf^] = 0, gives 



The last term can be turned into a surface integral over the boundary of the R'* (4D space-time region 
of integration). Since the initial and final field configurations are assumed given, 5$ = at the temporal 
beginning tin and end tout of this region, which implies that the surface term is zero. Factoring out the 5$ 
from the first two terms, and since the integral must vanish for arbitrary 5<^, we arrive at the Euler-lagrange 
equation of motion for a field. 



If the Lagrangian (density) C contains more fields, there is one such equation for each. The momentum density 
■k{x) of a field, conjugate to "I>(a::) is defined as: ^{x) = g-^^. 
For example, the standard electromagnetic action 

S = ~\ f (fx F^.F"", where F^, = d^A, - d,A^, 
4 Jk4 

gives the sourceless Maxwell's equations: 

d^F^- = 0, e^'""'9.F<,, = 0, 
where the field strength tensor F^i, and the Maxwell equations are invariant under the gauge transformations, 

A^ — > + d^€. 

The equations of motion of charged particles are given by the Lorentz-force equation, 

m = et Uv, 

dr 

where e is the charge of the particle and u'^ij) its four- velocity as a function of the proper time. 
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2. A complex-valued, quantum transition amplitude^ 

(Outi,,JK„,) := / P[$]e^^[*], (23) 



n 



^^The transition amplitude is closely related to partition function Z, which is a quantity that encodes the 
statistical properties of a system in thermodynamic equilibrium. It is a function of temperature and other 
parameters, such as the volume enclosing a gas. Other thermodynamic variables of the system, such as 
the total energy, free energy, entropy, and pressure, can be expressed in terms of the partition function or 
its derivatives. In particular, the partition function of a canonical ensemble is defined as a sum Z{l3) = 
e-'^^J, where /3 = l/(fcsr) is the ' inverse temperature', where T is an ordinary temperature and feg 
is the Boltzmann's constant. However, as the position and momentum pi variables of an ith particle in 
a system can vary continuously, the set of microstates is actually uncountable. In this case, some form of 
coarse-graining procedure must be carried out, which essentially amounts to treating two mechanical states 
as the same microstate if the differences in their position and momentum variables are 'small enough'. The 
partition function then takes the form of an integral. For instance, the partition function of a gas consisting 
of A'^ molecules is proportional to the 6 A'^— dimensional phase-space integral. 



/ d^p,d^a;'exph/3i?(p,;,x')], 
Jrbn 

where H = H{pi,x^), {i = 1, A^) is the classical Hamiltonian (total energy) function. 

Given a set of random variables Xi taking on values a;', and purely potential Hamiltonian function H{x^), 
the partition function is defined as 

Z(/3)=^exp [~pH{x')\. 

The function H is understood to be a real- valued function on the space of states {X\,X2,- ■ ■ } while /3 is a 
real- valued free parameter (conventionally, the inverse temperature). The sum over the x^ is understood to 
be a sum over all possible values that the random variable Xi may take. Thus, the sum is to be replaced by 
an integral when the Xi are continuous, rather than discrete. Thus, one writes 



Z{l3) = J dx'exp ^-l3H{x')^ , 



for the case of continuously- varying random variables Xi. 

Now, the number of variables Xi need not be countable, in which case the set of coordinates {x^} becomes 
a field (f) = (pix), so the sum is to be replaced by the Euclidean path integral (that is a Wick-rotated Feynman 
transition amplitude (|24p in imaginary time), as 

Z{^) = J V[^] exp [^H{^)]. 

More generally, in quantum field theory, instead of the field Hamiltonian H{(p) we have the action S{4>) of 
the theory. Both Euclidean path integral, 

Z{(f>) — J ©[(j!)] exp [—5(0)] , real path integral in imaginary time, (21) 

and Lorentzian one, 

Z((/)) = J T>[(l)] exp [15(0)] , complex path integral in real time, (22) 



are usually called 'partition functions'. While the Lorentzian path integral (|22p represents a quantum-field 
theory-generalization of the Schrodinger equation, the Euclidean path integral (I21|l represents a statistical- 
field-theory generalization of the Fokker-Planck equation. 
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where I' [4*] is 'an appropriate' Lebesgue-type measure, 

N 

AT— >oo ■ 



P[$] = Jim H^i, ii = l,...,n), 



so that we can 'safely integrate over a continuous spectrum and sum over a discrete spectrum 
of our problem domain il', of which the absolute square is the real-valued probability density 
function, 

P := \{Outt,Jlnt^J)f. 

This procedure can be redefined in a mathematically cleaner way if we Wick-rotate the 
time variable t to imaginary values, 1 1-^ t = t, thereby making all integrals real: 

|p[$]eiS[*] JEM^ |p[$]e-^[*]. (24) 

For example, in non-relativistic quantum mechanics, the propagation amplitude from Xa 
to Xb is given by the configuration path integral 

U{Xa,X,-T) = {x,\Xa) = {x,\ B-'"'^ \Xa) = J e'^I^^l , 

which satisfies the Schrodinger equation (in natural units) 



i-^U{xa,xi,;T) = HU{xa,Xb;T), where H = ---^ + V{xb). 



The phase-space path integral (without peculiar constants in the functional measure) 
reads 

rT 



U{qa,qb;T) = I n/ l^mmt)]]^^ {pd - H{q,p)) 



dt 



where the functions q{t) (space coordinates) are constrained at the endpoints, but the func- 
tions p{t) (canonically-conjugated momenta) are not. The functional measure is just the 
product of the standard integral over phase space at each point in time 



V[q{t)\V\p{t)\=\{^ [dq^dp, 



27r 

Applied to a non-relativistic real scalar field (j){x,t), this path integral becomes 
<0,(x,t)|e-^^|<^„(x,t)) = Jv[cj>]exp ij^ C{^)d'x , with Cicf>) = ^{d^cf>f - V icf>) . 
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2.2 Correlation functions and generating functional 

If we have two fields in the interacting theory, the corresponding two-point correlation func- 
tion, or two-point Green's function, is denoted by {n\T{(f){x)(f){y)}\Cl) , where the notation 
is introduced to denote the ground state of the interacting theory, which is generally 
different from |0) , the ground state of the free theory. The correlation function can be inter- 
preted physically as the amplitude for propagation of a particle or excitation between y and 
X. In the free theory, it is simply the Feynman propagator 

mw^mmoh^ = dh^ - „) = / ^^i^^^- 

We would like to know how this expression changes in the interacting theory. Once we have 
analyzed the two-point correlation functions, it will be easy to generalize our results to higher 
correlation functions in which more than two field operators appear. 
In general we have: 

(o\T{4>j{x)My)^wHI^TdtHim\o 

(n|r{<^(x)<^(y)}|Q) = hm \ 

T-oo(l-^e) (o\T{exp[-i J^^ dt Hi{t)]}\0 

where each diagram is built out of Feynman propagators, vertices and external points. 

The virtue of considering the time-ordered product is clear: It allows us to put everything 
inside one large T— operator. A similar formula holds for higher correlation functions of 
arbitrarily many fields; for each extra factor of (f) on the left, put an extra factor of (f)j on the 
right. 

In the interacting theory, the corresponding two-point correlation function is given by 

sum of all connected diagrams 
with two external points 



{i}\T{cj>ix)cj>{y)}\n) -- 
This is generalized to higher correlation functions as 

{n\T{^{xi)...4>{xn)m = 

In a scalar field theory, the generating functional of correlation functions is defined as 
Z[J] = j V[(t>] exp i j d^x [£ + J{x)(t>ix)] = {n\ e-'^^|0> = e-^-^'-^l; 



sum of all connected diagrams 
with n external points 
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this is a functional integral over (j){x) in which we have added a source term J{x)(j){x) to 
C = C{<t>). 

For example, the generating functional of the free Klein-Gordon theory is simply 



Z[J] = Zoexp 



^ j d^xct'yJ{x)DF{x-y)J{y) 



2.3 Quantization of the electromagnetic field 

Consider the path integral 

Z[A] = j VIA] e'-^I^l, where the action for the free e.-m. field is 

S[A] = j d^x -\{F^uf =\jd^xA^{x){d-'g''''-d>'d-)A,{x). 

Z[A] is the path integral over each of the four spacetime components: 

V[A] = V[AfV[A]'^V[AfV[Af. 

This functional integral is badly divergent, due to gauge invariance. Recall that F^^, and 
hence L, is invariant under a general gauge transformation of the form 

A^{x) Af,{x) + -dfj,a{x). 

e 

The troublesome modes arc those for which A^{x) = d^a[x), that is, those that are gauge- 
equivalent to A^j^{x) = 0. The path integral is badly defined because we are redundantly 
integrating over a continuous infinity of physically equivalent field configurations. To fix this 
problem, we would like to isolate the interesting part of the path integral, which counts each 
physical configuration only once. This can be accomplished using the Faddeev-Popov trick, 
which effectively adds a term to the system Lagrangian and after which we get 



exp 



d^x 



2C 



where ^ is any finite constant. 

This procedure needs also to be applied to the formula for the two-point correlation 
function 



{n\TO{A)\n) 



lim 

T^oo(l-ie) 



JV[A]0{A) exp 


i J^j, d'^x 


L 


/ P[A]exp 


i J^rp d^x £ 





which after Faddeev-Popov procedure becomes 



{?L\TO{A)\^) 



lim 











/I»[^]exp 


i J^j,d^x 


C 


-U^^A,)^' 
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3 Path-Integral TQFT 



3.1 Schwarz— type and Witten— type theories 

Consider a set of fields {0j} on a Riemannian n— manifold M (with a metric g^p) and real 
functional of these fields, S[(j)j\, which is the action of the theory. Also consider a set of 
operators Oa{(j)^) (labeled by some set of indices a), which are arbitrary functionals of the 
fields The vacuum expectation value (VEV) of a product of these operators is defined 

as the path integral (see [3B] ) 



A quantum field theory is considered topological if it possesses the following property: 



i.e., if the VEVs of some set of selected operators remain invariant under variations of the 
metric g^y on M. In this case, the operators Oa{4>i) called observables. 

There are two ways to formally guarantee that condition (j25p is satisfied. The first 
one corresponds to the situation in which both, the action, S, as well as the operators Oa, 
are metric independent. These TQFTs are called Schwarz-type. In the case of Schwarz- 
type theories one must first construct an action which is independent of the metric g^,^. The 
method is best illustrated by considering an example. Let us take into consideration the most 
interesting case of this type of theories: Chern-Simons gauge theory. The data in Chern- 
Simons gauge theory are the following: a differentiable compact 3-manifold M, a gauge group 
G, which will be taken simple and compact, and an integer parameter k. The action is the 
integral of the Chern-Simons form associated to a gauge connection A corresponding to the 
group G, 



Observables are constructed out of operators which do not contain the metric g^y 

In 

gauge invariant theories, as it is the case, one must also demand for these operators invariance 
under gauge transformations. An important set of observables in Chern-Simons gauge theory 
is constituted by the trace of the holonomy of the gauge connection A in some representation 





(25) 




(26) 
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R along a 1-cycle 7, that is the Wilson loopf^ 

Trji (Hol^(^)) = Tr^P exp / A. (27) 

The VEVs are labeled by representations Ri and embeddings 7^ of into M [55] 

(Tr/j^Pe^-^i ^ . . . TrR„Pe-^-^n " / t^^lTrfliPe-'^^i ^ • • • Tr^j^Pe-^-^n ^ei^^cs(A)_ 

A non-perturbative analysis of Chern-Simons gauge theory shows that the invariants associ- 
ated to the observables Oa{<Pi) are knot and link invariants of polynomial type as the Jones 
polynomial and its generalizations. The perturbative analysis has also led to this result and 
has shown to provide a very useful framework to study Vassiliev invariants. 

The second way to guarantee (j25p corresponds to the case in which there exists a sym- 
metry, whose infinitesimal form is denoted by (5, satisfying the following properties: 

(50„(,/.,) = 0, T^Mi) = &G^M^), (28) 
where T^,y(</>j) is the energy-momentum tensor of the theory, given by 

T^^Mi) = -^S[<^1 (29) 

while G^u{4>i) is some tensor. 

The fact that 5 in (j28p is a symmetry of the theory means that the transformations 5(1)^ 
of the fields are such that (^5'[</>j] = and 50a{4>i) — 0- Conditions (I28p lead formally to 
the following relation for VEVs: 

'X^ 0.2 ' ' ' ^cij,) — 



' P[<A,]O„,(0,)O,,(0J • • • Oa,(<Ajr^. exp (-5[0,]) (30) 
P[</)i]<5 {O^Mi)Oo.,{4>^) ■ ■ ■ Oa,(<^i)G^.exp {-S[(p,])) = 0, 



which implies that the quantum field theory can be regarded as topological. In (j30p it has 
been assumed that the action and the measure T>[(j)j\ are invariant under the symmetry 6. We 



^•^A holonomy on a smooth manifold is a general geometrical consequence of the curvature of the manifold 
connection, measuring the extent to which parallel transport around closed loops fails to preserve the geomet- 
rical data being transported. Related to holonomy is a Wilson loop, which is a gauge-invariant observable 
obtained from the holonomy of the gauge connection around a given loop. More precisely, a Wilson loop is 
a quantity defined by the trace of a path-ordered exponential of a gauge field transported along a closed 
curve (loop) 7, Wj = Tr(Pexp[i^ Af^dx'^]), where P is the path-ordering operator. 
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have assumed also in (j30p that the observables are metric-independent. This is a common 
situation in this type of theories, but it does not have to be necessarily so. In fact, in view 
of ([30]) . it would be possible to consider a wider class of operators satisfying: 

' = 50r(</'J, (31) 



where Oa'^((/)j) is a certain functional of the fields of the theory. 

This second type of TQFTs are called cohomological of Witten-type. One of its main rep- 
resentatives is Donaldson- Witten theory, which can be regarded as a certain twisted version 
of = 2 supersymmetric Yang-Mills theory. It is important to remark that the symmetry 6 
must be a scalar symmetry. The reason is that, being a global symmetry, the corresponding 
parameter must be covariantly constant and for arbitrary manifolds this property, if it is 
satisfied at all, implies strong restrictions unless the parameter is a scalar. 

Most of the TQFTs of cohomological type satisfy the relation: 5'[(/)j] = 6A{(j)^), for some 
functional A{(p^). This means that the topological observables of the theory (in particular the 
partition function itself) are independent of the value of the coupling constant. For example, 
consider the VEV ^ 

• • • O,^) = I V[ct>,]0^,{^,)0^,{<i>,) ■ ■ ■ O^^icPi) exp (-^5[<^.]) • (32) 

Under a change in the coupling constant, — A, one has (assuming that the 

observables do not depend on the coupling), up to first order in A: 



A I V[4>,]6 



Hence, observables can be computed either in the weak coupling limit, g ^ 0, or in the strong 
coupling limit, 5 — > oo. 



3.2 Hodge decomposition theorem 

The Hodge star operator * : ^^(M) — > il"~P(M), which maps any exterior p— form a E 
QP{M) into its dual (n — p)— form -ka £ f]"'~^(M) on a smooth n— manifold M, is defined as 
(see, e.g. [291 [30]) 

a A*/3 = /? A*a = (a,/3)^, * * a = (-l)P("-P)a, (for q, /3 G J1P(M)), 

The -k operator depends on the Riemannian metric g = gij on M and also on the orientation 
(reversing orientation will change the sign) [201 I21j . Using the star operator, for any two 
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p— forms a,(3 £ il^(M) with compact support on M we define bilinear and positive-definite 
Hodge L^— inner product as 

{a, (3) := / a A* p. (33) 

wliere q A * /5 is an n— form. 

Given the exterior derivative d : Qp(M) Q,p^^{M) on a smooth manifold M (see Ap- 
pendix), its Hodge dual (or, formal adjoint) is the codijferential 6, a linear map 6 : i}P{M) — > 
f)P~i(M), which is a generalization of the divergence, defined by [29\ I30j 

5 = (-l)"{p+i)+i*d^ so that d= (-!)"?'★ 5*. 

That is, if the dimension n of the manifold M is even, then 6 = — -k d-k. 
Applied to any p— form to G ^lf{M), the co differential 5 gives 

6u = (-l)"(P+^)+i*d*w, Jdc^ = {-l)''P+Udkdu;. 

If w = / is a 0— form, or function, then 6f = 0. If a p— form a is a codifferential of a 
(p + 1)— form /?, that is a = (5/3, then /? is called the coexact form. A p— form a is coclosed if 
(5a = 0; then *a is closed (i.e., d-ka = 0) and conversely. 

The Hodge codifferential 6 satisfies the following set of rules: 

• 5(5 = (5^ = 0, the same as dd = d"^ = 0; 

• 6k = {-1)P+^ k d; k6 = {-l)Pkd; 

• d6k: = kr 5d; kd6 = 6dk. 

The codifferential 6 can be coupled with the exterior derivative d to construct the Hodge 
Laplacian A : Qp{M) ilP(M), a harmonic generalization of the Laplace-Beltrami differ- 
ential operator, given by 

A = Sd + d5 = {d + 5f. 
A satisfies the following set of rules: 

(5 A = A (5 = 6d6; dA = Ad = d6d; kA = A-k . 

A p— form a is called harmonic iff 

Aa = <^=^ da = 5a = 0. 

Thus, a is harmonic in a compact domain D C Af iff it is both closed and coclosed in D. 
Informally, every harmonic form is both closed and coclosed. As a proof, we have: 

= (a, Aa) = (a, d6a) + (a, 6da) = {5a, 5a) + {da, da) . 
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Since /?) > for any form (3, {5a, 6a) and {da, da) must vanish separately. Thus, da = 
and 5a = 0. Ah harmonic p— forms on a smooth manifold M form the vector space H^{M). 

Now, the celebrated Hodge decomposition theorem (HDT) states that, on a compact ori- 
entable smooth n— manifold M (with n > p), any exterior form can be written as a unique 
sum of an exact form, a coexact form, and a harmonic form. More precisely, for any form 
uj £ J7^(M) there are unique forms a £ QP~^{M), (3 G Q,p^^{M) and a harmonic form 
7 G nP{M), such that 



HDT : 



any form 



exact 

da + 



coexact 



harmonic 

+ 7 



For the proof, see [291I30]. 

In physics community, the exact form da is called longitudinal, while the coexact form 5P is 
called transversal, so that they are mutually orthogonal. Thus any form can be orthogonally 
decomposed into a harmonic, a longitudinal and transversal form. For example, in fluid 
dynamics, any vector-field v can be decomposed into the sum of two vector-fields, one of 
which is divergence-free, and the other is curl-free. 

Since 7 is harmonic, dj = 0. Also, by Poincare lemma, d[da) = 0. In case w is a closed 
form, dio = 0, then the term 5(3 in HDT is absent, so we have the short Hodge decomposi- 
tion, 

UJ = da + j, (34) 

thus UJ and 7 differ by da. In topological terminology, uj and 7 belong to the same cohomology 
class [uj] G H'P{M). Now, by the de Rham theorems it follows that if C is any p— cycle, then 



Jc Jc 



c 

that is, 7 and uj have the same periods. More precisely, if uj is any closed form, then there 
exists a unique harmonic p— form 7 with the same periods as those of uj (see [29^ I31j). 

The Hodge-Weyl theorem |291 [30] states that every de Rham cohomology class has a 
unique harmonic representative. In other words, the space H^{M) of harmonic p— forms on 
a smooth manifold M is isomorphic to the pih de Rham cohomology group, 



ZP{M) _ Ker {d : nP{M) nP+^{M)) 
BPM ~ Im(d : ^P'^^M) nP{M)) 



(35) 



or, H-^{M) = H^j^{M). That is, the harmonic part 7 of HDT depends only on the global 
structure, i.e., the topology of M. 
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For example, in (2 + 1)D electrodynamics, jj— form Maxwell equations in the Fourier 
domain S are written as [32] 

dE = iujB, dB = 0, 

dH = -iuD + J, dD = Q, 

where H and uj are 0-forms (magnetizing field and field frequency), D (electric displacement 
field), J (electric current density) and E (electric field) are 1-forms, while B (magnetic field) 
and Q (electric charge density) are 2~forms. From d^ = it follows that the J and the Q 
satisfy the continuity equation 

dJ = iuQ. 

Constitutive equations, which include all metric information in this framework, are written 
in terms of Hodge star operators (that fix an isomorphism between p forms and (2 — p) forms 
in the (2 + 1) case) 

D = i.E, B = i.H. 
Applying HDT to the electric field intensity 1-form E, we get |33j 

E = d(l) + 5A + x, 

where (/; is a 0-form (a scalar field) and j4 is a 2-form; d(f) represents the static field and 5 A 
represents the dynamic field, and x represents the harmonic field component. If domain S is 
contractible, x is identically zero and we have the short Hodge decomposition, 

E = d(l) + 6 A. 

3.3 Hodge decomposition and gauge path integral 
3.3.1 Functional measure on the space of differential forms 

The Hodge inner product (j33p leads to a natural (metric-dependent) functional measure 
T)(j,[u}\ on VP[M), which normalizes the Gaussian functional integral 

yp^He'^'^l"^ = 1. (36) 

One can use the invariance of (|36p to determine how the functional measure transforms 
under the Hodge decomposition. Using HDT and its orthogonality with respect to the inner 
product (f33|) . it was shown in [35] that 

(w, u) = (7, 7) + {da, da) + {5(3, 5(3) = (7, 7) + (a, 5da) + {(3, d5(3) , (37) 

where the following differential/conferential identities were used |34j 

{da, da) = {a, 5 da) and {5(3,5(3) = {(3, d5(3). 
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Since, for any linear operator O, one has 

J Vfi[uj]expi{uj\Ou;) = det~^/^{0), 

(I36|) and §3) imply that 

3.3.2 Abelian Chern Simons theory 

Recall that the classical action for an Abelian Chern-Simons theory, 

S = [ AAdA, 
Jm 

is invariant (up to a total divergence) under the gauge transformation: 

A> — > A + df. (38) 
We wish to compute the partition function for the theory 

Z:= I^Vf.[AW^l^\ 

where Vq denotes the volume of the group of gauge transformations in (j38p , which must be 
factored out of the partition function in order to guarantee that the integration is performed 
only over physically distinct gauge fields. We can handle this by using the Hodge decompo- 
sition to parametrize the potential A in terms of its gauge invariant, and gauge dependent 
parts, so that the volume of the group of gauge transformations can be explicitly factored 
out, leaving a functional integral over gauge invariant modes only |35j . 
We now transform the integration variables: 

A I — > a, (3, J, 

where a, /3,7 parameterize respectively the exact, coexact, and harmonic parts of the con- 
nection A. Using the Jacobian (|37p as well as the following identity on 0-forms A = 5d, we 
get [35] 

Z = [ ^Vfi[a]Vn[P]Vfx[j] deti/2 (A) det^/^ (^j) ^iS ^ 
from which it follows that 

Vg= [ Vfi[a], (39) 
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while the classical action functional becomes, after integrating by parts, using the harmonic 
properties of 7 and the nilpotency of the exterior derivative operators, and dropping surface 
terms: 

S = -{p,i.6d6P) . 

Note that S depends only the coexact (transverse) part of A. Using ()39p and integrating over 
P yields: 

Z = J P^[7]det-i/2 (^^^5) _ 
Also, it was proven in [35] that 

det(*(5d(5) = det^/^ {{d6d) (6 d6)) = deti(d(5). 
As a consequence of Hodge duality we have the identity 

det (Sd) = det{d6), 

from which it follows that 

Z = I Vfi[j] det-3/4 (^Af,)) deti/2 (A) det^^ (^Af^)) . 

The operator A^^ is the transverse part of the Hodge Laplacian acting on 1— forms: 

Applying identity for the Hodge Laplacian A(p) [35] 

det (A(p)) = det det , 

we get 

det = det (A(i)) /det (A) 

and hence 



Z = J Vfi[-f] det-^/^ (A(i)) det^/^ (A) . 



The space of harmonic forms 7 (of any order) is a finite set. Hence, the integration over 
harmonic forms (|3.3.2p is a simple sum. 
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4 Non-Abelian Gauge Theories 

4.1 Intro to non-Abelian theories 

QED is the simplest example of a gauge theory coupled to matter based in the Abelian 
gauge symmetry of local U(l) phase rotations. However, it is possible also to construct 
gauge theories based on non-Abelian groups. Actually, our knowledge of the strong and weak 
interactions is based on the use of such non-Abelian generalizations of QED. 

Let us consider a gauge group G (see Appendix) with generators T", (a = 1, . . . , dimG) 
satisfying the Lie algebra 

A gauge field taking values on the Lie algebra of Q can be introduced = ^J^T", which 
transforms under a gauge transformations as 

where g is the coupling constant. The associated field strength is defined as 

F-, = d^Ai-d,Al-gr'''^A';^Al. 

Notice that this definition of the F^^^ reduces to the one used in QED in the Abelian case 
when /"'"^ = 0. In general, however, unlike the case of QED the field strength is not gauge 
invariant. In terms of F^,/ = F^^T"' it transforms as 

Fi,u UF^,U-\ 

The coupling of matter to a non-Abelian gauge field is done by introducing again the 
covariant derivative. For a field $ — ^ in a representation of Q, the covariant derivative 
is given by 

D^<^ = d^^ - igA^^T^<^. 

With the help of this we can write a generic Lagrangian for a non-Abelian gauge field coupled 
to scalars (f) and spinors ijj as 

1 



jC = + i^ppiP + D^c^D'^cP - V [Mi(</,) + 175^2(0)] V - V{<P). 

4.2 Yang-Mills theory 

In non-Abelian gauge theories, gauge fields are matrix-valued functions of space-time. In 
SU(N) gauge theories they can be represented by the generators of the corresponding Lie 
algebra, i.e., gauge fields and their color components are related by 

A,{x)=Al{x)^, (40) 
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where the color sum runs over the N"^ — 1 generators. The generators are hermitian, traceless 
N X N matrices whose commutation relations are specified by the structure constants /"''^ 
[391 



■ rabc ^ 

■' 2 



The normalization is chosen as 



Tr 



A'^ A^ 



Most of our applications will be concerned with SU (2) gauge theories; in this case the gen- 
erators are the Pauli matrices, 



A" 



with structure constants / 



abc ,abc 



Covariant derivative, field strength tensor, and its color components are respectively defined 

by 

D^ = d^ + igAf,, (41) 



1 



(42) 



The definition of electric and magnetic fields in terms of the field strength tensor is the same 
as in electrodynamics 



E'^ ix) 



-F' 



Oia 



x) 



B'" (x) 



(43) 



The dimensions of gauge field and field strength in 4D space-time are 

[A] = r\ [F]=r\ 

and therefore in absence of a scale, ~ M^^^^, with arbitrary constants M^^. In 

general, the action associated with these fields exhibits infrared and ultraviolet logarithmic 
divergencies. In the following we will discuss 



• Yang-Mills Theories: Only gauge fields are present. The Yang-Mills Lagrangian is 

Cym = -\f^""Ku = {F^'F^u) = i(E2 - b2). (44) 
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Quantum Chromodynamics: QCD contains besides the gauge fields (gluons), fermion 
fields (quarks). Quarks are in the fundamental representation, i.e., in SU(2) they are 
represented by 2-component color spinors. The QCD Lagrangian is (flavor dependences 
suppressed) 

^QCD = ^YM+^m, = {^^ D - m) tp , (45) 

with the action of the covariant derivative on the quarks given by 
{D^^y = {d^6^^ + igA^^) {iJ = l...N). 

Georgi-Glashow Model: In the Georgi-Glashow model [1^ (non-Abelian Higgs model), 
the gluons are coupled to a scalar, self-interacting (y((/>)) (Higgs) field in the adjoint 
representation. The Higgs field has the same representation in terms of the generators 
as the gauge field (jlO]) and can be thought of as a 3-component color vector in SU{2). 
Lagrangian and action of the covariant derivative are respectively 



Cgg = J^YM + C^, Cm = -D^cPD^cj) - VicP), (46) 



{D^<Pr = [D>.,^] " = {d^S^' - gr'^A'^^' . (47) 

4.2.1 Yang Mills action 

The general principle of least action, 

5S = 0, with S = J Cd^x, 

applied to the gauge fields, 

6Sym = -jd^'x Tr (F,..5F^") = " / Tr (^F^,^ [Z)^ 5^^) 
= 2 j d^xTT{6A''[D^',F^,]) 
gives the inhomogeneous field equations [39] 

[^M.i^n=r, (48) 

with j'^ the color current associated with the matter fields 
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For QCD and the Georgi-Glashow model, these currents are given respectively by 

As in electrodynamics, the homogeneous field equations for the Yang-Mills field strength 

0, 



with the dual field strength tensor 



are obtained as the Jacobi identities of the covariant derivative, 

[D^, [D,, D,]] + [D,, [Dp, D,,]] + [D„ [D,, D^]] = 0. 

4.2.2 Gauge transformations 

Gauge transformations change the color orientation of the matter fields locally, i.e., in a 
space-time dependent manner, and are defined as 



U (x) = exp {iga (x)} = exp < iga 



r"1 



with the arbitrary gauge function (x). Matter fields transform covariantly with U 

ipvUip, <p — >U(l)Ul (51) 

The transformation property of A is chosen such that the covariant derivatives of the matter 
fields D^ip and Df^(j) transform as the matter fields ip and (p respectively. As in electrody- 
namics, this requirement makes the gauge fields transform inhomogeneously [39] 



(x) U (x) (^A^ (x) + la^) [/t (x) = (a 



(52) 



resulting in a covariant transformation law for the field strength 

i^M- UF^M^. (53) 
Under infinitesimal gauge transformations {{ga"" (x) | <C 1) 

(x) Al (x) - d^a^ (x) - gr'^'a" (x) (x) . (54) 
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As in electrodynamics, gauge fields which are gauge transforms of = are called pure 
gauges and are, according to ([52]) . given by 

AP^x) = Uix)^d^uHx) ■ (55) 

Physical observables must be independent of the choice of gauge (coordinate system in color 
space). Local quantities such as the Yang-Mills action density Tr (F^'^(a;)F^jy(x)) or mat- 
ter field bilinears like 'ijj{x)ip{x),(l)"'{x)(j)'^{x) are gauge invariant, i.e., their value does not 
change under local gauge transformations. One also introduces non-local quantities which, 
in generalization of the transformation law (j53p for the field strength, change homogeneously 
under gauge transformations. In this construction a basic building block is the path ordered 
integral 

n{x,y,C) = Pexp!^-ig da^A^, {x{a))^ = P exp !^-ig j dx^'A^. (56) 

It describes a gauge string between the space-time points x = x{sq) and y = x{s). 17 satisfies 
the differential equation 

d9. . dx^ . ^ 

- = (57) 

Gauge transforming this differential equation yields the transformation property of Vt 

n{x,y,C) ^u{x)n{x,y,c)uHy)- (58) 

With the help of 0,, non-local, gauge invariant quantities like 

Tr {F^'''{x)n {x, y, C) F^,{y)) , ^{x)n (x, y, C) ^{y), 
or closed gauge strings, the following SU(N)-Wilson loops 

Wc = ^T,{n{x,x,C)) (59) 



can be constructed. For pure gauges (j55|) . the differential equation (j57|) is solved by 

nPnx,y,C) = U{x)U\y). (60) 

While ^p{x)^} {x,y,C) ip{y) is an operator which connects the vacuum with meson states for 
SU{2) and SU{3), fermionic baryons appear only in SU{3) in which gauge invariant states 
containing an odd number of fermions can be constructed. In SU(3) a point-like gauge 
invariant baryonic state is obtained by creating three quarks in a color antisymmetric state 
at the same space-time point 

V'(x) ~ e''^''^P''{x)iJ\x)i;''{x). 
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Under gauge transformations, 



V^(x) e''''Uaa{xW{x)U,,p{xW\x)U,^{x)i;\x) 
= det {U{x)) e''^''i:''{x)il:^{x)^'\x) . 

Operators that create finite size baryonic states must contain appropriate gauge strings as 
given by the fohowing expression 

^(x, y, z) ~ e"^^[r!(^/, x, Ci)^(x)]" [^(n, y, C2)^{y)f [n{u, z, C^)^{z)r • 

The presence of these gauge strings makes ij) gauge invariant as is easily verified with the help 
of the transformation property (j58p . Thus, gauge invariance is enforced by color exchange 
processes taking place between the quarks. 

4.3 Quantization of Yang— Mills theory 

Gauge theories are formulated in terms of redundant variables. Only in this way, a covariant, 
local representation of the dynamics of gauge degrees of freedom is possible. For quantiza- 
tion of the theory both canonically or in the path integral, redundant variables have to be 
eliminated. This procedure is called gauge fixing. It is not unique and the implications of a 
particular choice are generally not well understood. In the path integral one performs a sum 
over all field configurations. In gauge theories this procedure has to be modified by making 
use of the decomposition of the space of gauge fields into equivalence classes, the gauge or- 
bits. Instead of summing in the path integral over formally different but physically equivalent 
fields, the integration is performed over the equivalence classes of such fields, i.e., over the 
corresponding gauge orbits. The value of the action is gauge invariant, i.e., the same for all 
members of a given gauge orbit. Therefore, the action is seen to be a functional defined on 
classes (gauge orbits) [39j. Also the integration measure 

V[A]= X{dAl{x). 

x,fi,a 

is gauge invariant since shifts and rotations of an integration variable do not change the value 
of an integral. Therefore, in the naive path integral 

Z[A] = J P[^]e'^[^] oc JlldU{x). 

a 'volume' associated with the gauge transformations dU (x) can be factorized and 
thereby the integration be performed over the gauge orbits. To turn this property into 
a working algorithm, redundant variables are eliminated by imposing a gauge condition, 
f[A] = 0, which is supposed to eliminate all gauge copies of a certain field configuration A. 
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In other words, the functional / has to be chosen such that, for arbitrary field configurations, 
the equation, /[Al^l] = 0, determines uniquely the gauge transformation U. If successful, 
the set of all gauge equivalent fields, the gauge orbit, is represented by exactly one represen- 
tative. In order to write down an integral over gauge orbits, we insert into the integral the 
gauge-fixing (5— functional 

S[f{A)] = ll n Hr{A{x))]. 

X a=l 

This modification of the integral however changes the value depending on the representative 
chosen, as the following elementary identity shows 

6 {x — a) 



W (a) 



g{a) = 0. 



This difficulty is circumvented with the help of the Faddeev-Popov determinant A / [A] 
defined implicitly by 

A;[^] I P[C/]<5[/(^[^l)] =1. 

Multiplication of the path integral Z[A] with the above "1" and taking into account the gauge 
invariance of the various factors yields 

Z[A] = j V[U] j V [A] e'-^l^lA/ [A] 5 [/ (^ylt^l) 

= y P [[/] y P [A\ e^^I^''"'] A^ [^[^]] 5 [f (^[^])" . 

The gauge volume has been factorized and, being independent of the dynamics, can be 
dropped. In summary, the final definition of the generating functional for gauge theories is 
given in terms of a sum over gauge orbits, 

Z[J]= Jv [A] Af [A] S (/ [A] ) e'^H+'/<i'^-^''^M. 



4.3.1 Faddeev— Popov determinant 

For the calculation of A/ [A] , we first consider the change of the gauge condition [^4] under 
infinitesimal gauge transformations. Taylor expansion 



Jx 



Aim 



rxiAHjd'yJl^^Kiy) 

[A] + j d^yY,M{x,y- a, b) a' (y) , 



49 



with given by infinitesimal gauge transformations, 



Alix) ^ Al{x)-d,a-{x)-gr''^a'ix)Al{x) 



yields 



M {x,y;a,b) 



d,5'^' + gr^AUy) 



{yy 



In the second step, we compute the integral 

by expressing the integration P \U] as an integration over the gauge functions a. We finally 
change to the variables (5 = Ma, 

[A] = I det M\-^ Jv [f3] 6 [f {A) - f3] , 

and arrive at the final expression for the Faddeev-Popov determinant [39] 

Af [A] = |detM|. 

Examples: 

• Lorentz gauge 



fS (A) 
M{x,y;a,b) 



d^^A';, (x)-x" (x), 

's-ba_gfabc^c (y)0A 5(4) 



• Coulomb gauge 



M{x,y;a,h) 



divA'^ (x) 
5'^'A + gr'''A^y)Vy)5(^^ {x - y) 



Axial gauge 



f^{A) = n^^Al{x)-x''{x), 
M{x,y;a,b) = -6''^n^d^5^'^'^ {x - y) 
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4.4 Basics of Conformal field tlieory 

A conformal field theory (CFT) is a quantum field theory (or, a statistical mechanics model 
at the critical point) that is invariant under conformal transformations. Conformal field 
theory is often studied in 2D where there is an infinite-dimensional group of local conformal 
transformations, described by the holomorphic functions. CFT has important applications in 
string theory, statistical mechanics, and condensed matter physics. For a good introduction 
to CFT see dH US]. We consider here only chiral CFTs in 2D (see [i3|), where 'chiral0 
means that all of our fields will be functions of a complex number z = x + iy only and not 
functions of its conjugate z. 

To formally describe a 2D CFT we give its 'conformal data', including a set of primary 
fields, each with a conformal dimension A, a table of fusion rules of these fields and a central 
charge c. Data for three CFTs are given in Table [H 

The operator product expansion (OPE) describes what happens to two fields when their 
positions approach each other. We write the OPE for two arbitrary fields (/>j and <j)j as 

lim (/),(z)0,(u;) = C^^{z - «;)^'=-^-^.- ^,{w), 

k 

where the structure constants C^j are only nonzero as indicated by the fusion table. Note 
that the OPE works inside a correlator. For example, in the Z3 para-fermion CFT (see Table 
[T]), since cti x tpi = e, for arbitrary fields (j)^ we have [43] 



lim {(l)^{zi) . . . (t)j^i{zM) <Ti{z)i:^{w)) ~ 

2:— >ui 

_ u;)2/5-i/i5-2/3^^^(^^) ^,,{zM)e{w)). 

In addition to the OPE, there is also an important 'neutrality' condition: a correlator 
is zero unless all of the fields can fuse together to form the identity field 1. For example, 
in the Z3 para-fermion field theory (V'2V'i) / since V'2 ^ V'l = 1; but {ipiipi) = since 

V'l X ^1 = V'2 / 1- 

Let us look at what happens when a fusion has more than one possible result. For example, 
in the Ising CFT, a x a = 1 + ip. Using the OPE, we have 

lim a{wi)a{w2)-^ -. — -rj^+ {wi - W2f^^ ip, (61) 

wi^W2 [Wi — W2) ' 

where we have neglected the constants C^j. If we consider (era), the neutrality condition picks 
out only the first term in (j6ip where the two cr's fuse to form 1. Similarly, (aaip) results in 
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In general, a chiral field is a holomorphic field W(z) which transforms as 

L„W(z) = - (n + l)A2"W(z), with Lr,W(z) = Q, 
and similarly for an anti-chiral field. Here, A is the conformal weight of the chiral field W . 
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Ising CFT: (c= 1/2) 



Chiral Bose Vertex: (c = 1) 
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Z3 Parafermion CFT: 
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0-1 


e 


£72 


(^2 + V^i 






<72 


CTl 


e 


1 + e 


CJl + V2 
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0-2 
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CJi + V'2 


0-2 + V'l 


l + € 



Table 1: Conformal data for three CFTs. Given is the hst of primary fields in the CFT with 
their conformal dimension A, as weh as the fusion table. In addition, every CFT has an 
identity field 1 with dimension A = which fuses trivially with any field (1 x 0^ = for any 
(?!)j). Note that fusion tables are symmetric so only the lower part is given. In the Ising CFT 
the field tp is frequently notated as e. This fusion table indicates the nonzero elements of the 
fusion matrix N^j^. For example in the Z3 CFT, since (Ji x (72 = 1 + e, iV^jo-j ~ -^^10-2 ~ 
and -^CTi<T2 ~ ^ ^ equal to 1 or e. 
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the second term of (j6ip where the two cj's fuse to form ^ which then fuses with the additional 
-0 to make 1. 

Fields may also fuse to form the identity in more than one way. For example, in the 
correlator {a{wi)a{w2)cr{w3)a{w4)) of the Ising CFT, the identity is obtained via two possible 
fusion paths — resulting in two different so-called 'conformal blocks'. On the one hand, one 
can fuse cr^wi) and (t('W2) to form 1 and similarly fuse (7(103) (^{104^) to form 1. Alternately, 
one can fuse cr{wi) and (7(102) to form ip and fuse (7(103) and (t(w4) to form ip then fuse the 
two resulting ^ fields together to form 1. The correlator generally gives a linear combination 
of the possible resulting conformal blocks. We should thus think of such a correlator as living 
in a vector space rather than having a single value. (If we instead choose to fuse 1 with 3, 
and 2 with 4, we would obtain two blocks which are linear combinations of the ones found 
by fusing 1 with 2 and 3 with 4. The resulting vectors space, however, is independent of 
the order of fusion). Crucially, transporting the coordinates Wi around each other makes a 
rotation within this vector space. 

To be clear about the notion of conformal blocks, let us look at the explicit form of the 
Ising CFT correlator [53] 

lim {a(0)a(z)a(l)a(w)) = a+ F+ + a_ F_, 

ui— >oo 

F±(z) ~ (wz(l - z))-'^/^^Jl±^/l- z, 

where a+ and a_ are arbitrary coefficients. When z — > we have ~ z~^^^ whereas 
F- ~ z^^^. Comparing to (j6ip we conclude that is the result of fusing (7(0) x a(z) 1 
whereas F_ is the result of fusing (7(0) x (7(z) ip. As z is taken in a clockwise circle 
around the point z = 1, the inner square-root changes sign, switching and Thus, this 
'braiding' (or 'monodromy') operation transforms 

Having a multiple valued correlator (I.e., multiple conformal blocks) is a result of having 
such branch cuts. Braiding the coordinates (w's) around each other results in the correlator 
changing values within its allowable vector space. 

A useful technique for counting conformal blocks is the Bratteli diagram. In Figure [T] we 
give the Bratteli diagram for the fusion of multiple cr fields in the Ising CFT. Starting with 
1 at the lower left, at each step moving from the left to the right, we fuse with one more a 
field. At the first step, the arrow points from 1 to cr since 1 x cr = o". At the next step a 
fuses with cr to produce either ^ or 1 and so forth. Each conformal block is associated with 
a path through the diagram. Thus to determine the number of blocks in {aaaa) we count 
the number of paths of four steps in the diagram starting at the lower left and ending at 1. 

A particularly important CFT is obtained from a free Bose field theory in (1+1)D by 
keeping only the left moving modes. The free chiral Bose field (t)(z), which is a sum of left 
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tj a fj . . . 

1 V 1 

Figure 1: Bratteli diagram for fusion of multiple cJi fields in the Z3 para-fermion CFT (mod- 
ified and adapted from [i3]). 



moving creation and annihilation operators, has a correlator {(p{z)(f){z')) = — log{z — z'). We 
then define the normal ordered 'chiral vertex operator', e^"'^^^^ which is a conformal field. 
Since </> is a free field. Wick's theorem can be used to obtain 



4.5 Chern— Simons theory and Jones polynomial 
4.5.1 Chern Simons theory and link invariants 

Consider SUiT)]^ non-Abelian Chern-Simons theory 



Scs\a\ = T- /" tr ( a A da + I 



2 > 

-a A a A a 



We modify the action by the addition of sources, j'^-, according to £ ^ £-|-tr (j • a). We take 
the sources to be a set of particles on prescribed classical trajectories. The particle carries 
the spin jj representation of SU{2). As we have seen above, there are only A: + 1 allowed 
representations; later in this subsection, we will see that if we give a particle a higher spin 
representation than j = k/2, then the amplitude will vanish identically. Therefore, ji must 
be in allowed set of k + 1 possibilities: 0, ^, . . . , ^. The functional integral in the presence of 
these sources can be written in terms of Wilson loops, 1^7j,jJo], which are defined as follows. 
The holonomy U^j [a] is an SU (2) matrix associated with a curve 7. It is defined as the 
path-ordered exponential integral of the gauge field along the path 7 [43j 



'2 • • • 



C/^,[a] = Pe'^^"^^^-''' = Vi"/ dsi ds2 

n=0 

/ dsn [7(si) • (7(51)) T^'... iisn) ■ a^" (7(sn)) T^-] , 
Jo 

where V is the path-ordering symbol. The Lie algebra generators T- are taken in the spin j 
representation. j{s), s £ [0, 27r] is a parametrization of 7; the holonomy is clearly independent 
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of the parametrization. The Wilson loop is the trace of the holonomy: 

W^j[a]=ti{U^j[a]). (62) 

Let us consider the simplest case, in which the source is a quasi-particle-quasi-hole pair of 
type j which is created out of the ground state, propagated for a period of time, and then 
annihilated, returning the system to the ground state. The transition amplitude for such a 
process is given by the path integral: 

(0|0),,,. = I P[a]ei^^^Ht^^^^.[a], 

where 7 is the spacetime loop formed by the trajectory of the quasi-particle-quasi-hole pair. 
The Wilson loop was introduced as an order parameter for confinement in a gauge theory 
because this amplitude roughly measures the force between the quasi-particle and the quasi- 
hole. If they were to interact with a confining force V{r) ~ r, then the logarithm of this 
amplitude would be proportional to the the area of the loop; if they were to have a short- 
ranged interaction, it would be proportional to the perimeter of the loop. However, Chern- 
Simons theory is independent of a metric, so the amplitude cannot depend on any length 
scales. It must simply be a constant. For j = 1/2, we will call this constant d. As the 
notation implies, it is, in fact, the quantum dimension of a j = 1/2 particle. As we will see 
below, d can be determined in terms of the level k, and the quantum dimensions of higher 
spin particles can be expressed in terms of d. 

We can also consider the amplitude for two pairs of quasi-particles to be created out of 
the ground state, propagated for some time, and then annihilated, returning the system to 
the ground state [43] 

(010)71 Ji;7„i2 = / ^[«]e'^^«t'^lW^^,,HM^y,,v[a]. 

This amplitude can take different values depending on how 7 and 7' are linked. If the curves 
are unlinked the integral must give d^, but when they are linked the value can be nontrivial. 
In a similar way, we can formulate the amplitudes for an arbitrary number of sources. 

It is useful to think about this history as a two step process: from t = —00 to t = and 
from t = to t = 00 (the two pairs are created at some time t < and annihilated at some 
time t > 0). At t = 0~, the system is in a four-quasi-particle stateEl Let us call this state i/j: 

m] = ! ^[a(x,t)] : W^, j[a] : T^y ,,-,[a] e^--''*^'''"^^cs^ 

./a{x,0)=A{x) 

^^Quasiparticles and quasi-holes are topologically equivalent il G — SU{2), so we will use 'quasi-particle' to 
refer to both. 
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where 7_ and are the arcs given by ^{t) and 7'(t) for t < 0. ^(x) is the value of the 
gauge field on the t = spatial slice; the wave- functional il^[A\ assigns an amplitude to every 
spatial gauge field configuration. For G = SU{2) and k > 1, there are actually two different 
four-quasi-particle states: if particles 1 and 2 fuse to the identity field j = 0, then particles 
3 and 4 must as well; if particles 1 and 2 fuse to j = 1, then particles 3 and 4 must as 
well. These are the only possibilities0 Which one the system is in depends on how the 
trajectories of the four quasi-particles are intertwined. Although quasi-particles 1 and 2 were 
created as a pair from the vacuum, quasi-particle 2 braided with quasi-particle 3, so 1 and 2 
may no longer fuse to the vacuum. In just a moment, we will see an example of a different 
four-quasi-particle state. 

We now interpret the t = to t = oo history as the conjugate ofai = — (X)tot = 
history. In other words, it gives us a four quasi-particle bra rather than a four quasi-particle 
ket [43j 

X*[A] = [ P[a(x,t)] : W,^M : W^y ^,.,[a] e/o°='^*/'^'-^cs. 

Ja(x,0)=A(x) + 

In the state |x), quasi-particles 1 and 2 fuse to form the trivial quasi-particle, as do quasi- 
particles 3 and 4. Then we can interpret the functional integral from t = — oo to t = cxd as 
the matrix element between the bra and the ket: 

(X|V') = J ^[«]e'^^«[»]vr^,jJa]W^„,,[a]. 

Now, observe that \ip) is obtained from \x) by taking quasi-particle 2 around quasi-particle 
3, i.e., by exchanging quasi-particles 2 and 3 twice, {i/j) = ^(cl) |x)- Hence, 

(xIpH) \x) = I V[a] e'^^^W W.^jAa] W^^.M- (63) 

In this way, we can compute the entries of the braiding matrices p(o'j) by computing functional 
integrals such as the one on the right-hand-side of (I63p . 

Consider, now, the state p{cr2) \x)j iii which particles 2 and 3 are exchanged just once. 
From the figure, we see that 

{X\p{cr2)\x) = d, {X\p{cr2^) \x) = d, 

since both histories contain just a single unknotted loop. Meanwhile, (xlx) = d'^. 

Since the four-quasi-particle Hilbert space is 2D, p{a2) has two eigenvalues, Ai, A2, so 
that 

p{a) - (Ai + A2) + Xi\2p{(J-^) = 0. 
^^For k = 1, fusion to j = 1 is not possible. 
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Taking the expectation value in the state we find: 

a — (Ai + A2) a + AiA2a = 0, so that d, = — — . 

Ai + A2 

Since the braiding matrix is unitary, Ai and A2 are phases. The overall phase is unimportant 
for quantum computation, so we really need only a single number. In fact, this number can 
be obtained from self-consistency conditions [H]. However, the details of the computation of 
Ai, A2 within is technical and requires a careful discussion of framing; the result is [i2\ that 

Ai = _e-3W2(fc+2)^ ^ g7ri/2(fc+2)^ 
11 11 

These eigenvalues are simply i^Q ' ^ = Xi, Rf'^ = A2. Consequently, 

d = 2 cos (^Y^^ and (64) 

q-'^''p{c7,)-q'/^p{ar')=q-q~\ (65) 

where q = . Since this operator equation applies regardless of the state to which it 

is applied, we can apply it locally to any given part of a knot diagram to relate the amplitude 
to the amplitude for topologically simpler processes. This is an example of a skein relation; in 
this case, it is the skein relation which defines the Jones polynomial. In arriving at this skein 
relation, we are retracing the connection between Wilson loops in Chern-Simons theory and 
knot invariants which was made in the remarkable paper [12]. In this paper, Witten showed 
that correlation functions of Wilson loop operators in SU{2)k Chern-Simons theory are equal 
to corresponding evaluations of the Jones polynomials^ which is a topological invariant of 
knot i/ieorjll 

j V[a] VF^^^i [a]... W^^i_ [a] e^^^^W = Viiq). (66) 

^^A link is a finite family of disjoint, smootli, oriented or un-oriented, closed curves in or equivalently 
S^. A knot is a link with one component. The Jones polynomial Vhit) is a Laurent polynomial in the variable 
\/t, which is defined for every oriented link L but depends on that link only up to orientation preserving 
diffeomorphism, or equivalently isotopy, of R^. Links can be represented by diagrams in the plane. The Jones 
polynomial of a knot (and generally a link with an odd number of components) is a Laurent polynomial in t 

m- 

^*Knot theory is the area of topology that studies mathematical knots. Formally, a knot is an embedding 
of a circle in 3D Euclidean space, R"^. Two knots are equivalent if one can be transformed into the other 
via a deformation of R^ upon itself (known as an ambient isotopy). These transformations correspond to 
manipulations of a knotted string that do not involve cutting the string or passing the string through itself. 
Knots can be described in various ways. Given a method of description, however, there may be more than 
one description that represents the same knot. For example, a common method of describing a knot is a 
planar knot diagram. However, any given knot can be drawn in many different ways using a knot diagram. 
Therefore, a fundamental problem in knot theory is determining when two descriptions represent the same 
knot. One way of distinguishing knots is by using a knot invariant, a quantity which remains the same even 
with different descriptions of a knot. The concept of a knot has also been extended to higher dimensions by 
considering mD spheres S"" in R". 
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VL{q) is the Jones polynomial associated with the link L = U . . . U 7„, evaluated at 
q = using the skein relation (j65p . Note that we assume here that all of the quasi- 

particles transform under the j = \ representation of SU{2). The other quasi-particle types 
can be obtained through the fusion of several j = 1/2 quasi-particles. 

5 Appendix 

5.1 Manifolds and bundles 

Geometrically, a manifold is a nonlinear (i.e., curved) space which is locally homeomorphic 
(i.e., topologically equivalent) to a linear (i.e., flat) Euclidean space M"; e.g., in a magnifying 
glass, each local patch of the apple surface looks like a plane, although globally (as a whole) 
the apple surface is totally different from the plane. Physically, a configuration manifold is a 
set of all degrees of freedom of a dynamical system. 

More precisely, consider Consider a set M (see Figure [2]) which is a candidate for a 
manifold. Any point x G iW0 has its EucMean chart; given by a 1-1 and onto map (p^ : 
M — > M", with its Euclidean image Vi = (pi{Ui). More precisely, a chart ip^ is defined by 

(fi:MDU^3x^ ip^{x) eViC M", 

where Ui <Z M and Vi C are open sets. 




Figure 2: Geometric picture of the manifold concept. 

^^Note that sometimes we will denote the point in a manifold M by m, and sometimes by x (thus implicitly 
assuming the existence of coordinates x = (a;'))- 
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Clearly, any point x £ M can have several different charts (see Figure [2|). Consider a case 
of two charts, ip^, (pj : M ^ M", having in their images two open sets, Vij = ^Pj_{Ui n Uj) and 
Vji = ^j{Ui nUj). Then we have transition functions ip^j between them. 



-1 



locally given by (x) = (/j^ (93 • ^ (x) ) . 



If transition functions (/^j^ exist, then we say that two charts, ^p^ and ^pj are compatible. 
Transition functions represent a general (nonlinear) transformations of coordinates, which 
are the core of classical tensor calculus. 

A set of compatible charts ip^ : M ^ M", such that each point x G M has its Euclidean 
image in at least one chart, is called an atlas. Two atlases are equivalent iff all their charts 
are compatible (i.e., transition functions exist between them), so their union is also an atlas. 
A manifold structure is a class of equivalent atlases. 

Finally, as charts (p^ : M ^ were supposed to be 1-1 and onto maps, they can be either 
homeomorphisms, in which case we have a topological (C") manifold, or diffeomorphisms, in 
which case we have a smooth (C^) manifold. 

On the other hand, tangent and cotangent bundles, TM and T*M, respectively, of a 
smooth manifold M, are special cases of a more general geometrical object called fibre bundle, 
where the word fiber y of a map tt : Y ^ X denotes the preimage 7r~^(x) of an element 
X G X. It is a space which locally looks like a product of two spaces (similarly as a manifold 
locally looks like Euclidean space), but may possess a different global structure. To get a 
visual intuition behind this fundamental geometrical concept, we can say that a fibre bundle 
y is a homeomorphic generalization of a product space X xV (see Figure [3|), where X and V 
are called the base and the fibre, respectively, vr : y — > X is called the projection, = tt^^{x) 
denotes a fibre over a point x of the base X, while the map / = vr^^ : X ^ Y defines the 
cross-section, producing the graph (x, /(x)) in the bundle Y (e.g., in case of a tangent bundle, 
/ = X represents a velocity vector-field, so that the graph in a the bundle Y reads (x,x)). 



bundle Y'^XxV 




V 



Y 








f\ 1 1^ 



X 



X 



Figure 3: A sketch of a fibre bundle y ~ X x as a generalization of a product space X xV; 
left - main components; right ~ a few details (see text for explanation). 



59 



A principal G— bundle is a bundle n : Y — > X generated by a Lie group G (see below) 
such that the group G preserves the fibers of the bundle Y. 

The main reason why we need to study fibre bundles is that all dynamical objects (in- 
cluding vectors, tensors, differential forms and gauge potentials) are their cross-sections, 
representing generalizations of graphs of continuous functions. For more technical details, 
see [201 m]. 

5.2 Lie groups 

A Lie group is a both a group and a manifold. More precisely, a Lie group is a smooth 
manifold M that has at the same time a group G— structure consistent with its manifold 
M— structure in the sense that group multiplication /j, : G x G ^ G, {g,h) gh and 
thegroup inversion u : G ^ G, g ^ g^^ are smooth functions. A point e S G is called the 
group identity element. For any Lie group G in a neighborhood of its identity element e it 
can be expressed in terms of a set of generators (a = 1, . . . , dim G) as 

D{g) = exp[-ia„r"] = V ^^a^^ . . . a.^T'^^ . . . T"", 

n=0 

where € C are a set of coordinates of M in a neighborhood of e. Because of the general 
Baker-Campbell-Haussdorf formula, the multiplication of two group elements is encoded in 
the value of the commutator of two generators, that in general has the form 

l^^a i jabcrpc 

where f"'^'^ G C are called the structure constants. The set of generators with the commutator 
operation form the Lie algebra associated with the Lie group. Hence, given a representation of 
the Lie algebra of generators we can construct a representation of the group by exponentiation 
(at least locally near the identity). 

In particular, for 5C/(2)— group, each group element is labeled by three real numbers a^, 
{k = 1,2,3). We have two basic representations: one is the fundamental representation (or 
spin i) defined by 

Z;i(a,,) = e-W-', 

2 

with o"' the Pauli matrices. The second one is the adjoint (or spin 1) representation which 
can be written as 

Di{ak) = e''""-^" , where 

/000\ /00-l\ /010\ 

J^= 1, j2= 00 , J^= -100 . 

Vo-10/ Vloo/ \ooo/ 
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Actually, J generate rotations around the x, y and z axis respectively. 

Let M be a smooth manifold. An action of a Lie group G (with the unit element e) on 
M is a smooth map <j) : G x M ^ M, such that for all x G M and g,h € G, (i) <j){e, x) = x 
and (ii) (j) {g,(j){h,x)) = (j){gh,x). For more technical details, see |2m [2T]. 

5.3 Differential forms and Stokes theorem 

Given the space of exterior differential p— forms Qp(M) on a smooth manifold M, we have 
the exterior derivative operator d : Q{M) — > Q,p~^^{M) which generalizes ordinary vector 
differential operators {grad, div and curl see |29[ [5T| l2Uj ) and transforms p— forms lo into 
(p + 1)— forms did, with the main property: dd = d^ = 0. Given a p— form a G ^^{M) and a 

form (3 G r2'^(M), their exterior product is a (p + g)— form a f\ (3 ^ 0^+^(M), where A is 
their anti-commutative exterior (or, 'wedge') product. 

As differential forms are meant for integration, we have a generalization of all integral 
theorems from vector calculus in the form of the Stokes theorem: for the p— form in an 
oriented nD domain C, which is a p— chain with a (p — 1)— boundary dC, 



For any p— chain on a manifold M, the boundary of a boundary is zero [18j, that is, ddC = 



A p— form P is called closed if its exterior derivative d = didx^ is equal to zero, dP = 
O.From this condition one can see that the closed form (the kernel of the exterior deriva- 
tive operator d) is conserved quantity. Therefore, closed p— forms possess certain invariant 
properties, physically corresponding to the conservation laws (see e.g., [36} [2T]). 

Also, a p— form /? that is an exterior derivative of some (p — 1)— form a, f3 = da, is called 
exact (the image of the exterior derivative operator d). By Poincare lemma, exact forms 
prove to be closed automatically, dp = d{da) = 0. 

bmce every exact form is closed. The converse is only partially true, by Poincare 

lemma: every closed form is locally exact. In particular, there is a Poincare lemma for 
contractible manifolds: Any closed form on a smoothly contractible manifold is exact. The 
Poincare lemma is a generalization and unification of two well-known facts in vector calculus: 
(i) If curli^ = 0, then locally F = grad/; and (ii) If divF = 0, then locally F = curlG. 

A cycle is a p— chain, (or, an oriented p— domain) C G Cp{M) such that dC = 0. A 
boundary is a chain G such that C = dB, for any other chain B G Cp{M). Similarly, a 
cocycle (i.e., a closed form) is a cochain uj such that duj = 0. A coboundary (i.e., an exact 
form) is a cochain lo such that u = d9, for any other cochain 9. All exact forms are closed 
(lo = d9 ^ duj = 0) and all boundaries are cycles (C = dB ^ dC = 0). Converse is true 
only for smooth contractible manifolds, by Poincare lemma. 




(67) 



^2 = 0. 
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Integration on a smooth manifold M should be thought of as a nondegenerate bilinear 
pairing (,) between p— forms and p— chains (spanning a finite domain on M). Duality of 
p— forms and p— chains on M is based on the de Rham's 'period', defined as [29^ I34j 



Period := 




:=(C7,^), 



where C is a cycle, w is a cocycle, while {C,uj) = uj{C) is their inner product (C, w) : 
QP{M) X Cp{M) — > M. From the Poincare lemma, a closed p— form cu is exact iff {C,uj) = 0. 

The fundamental topological duality is based on the Stokes theorem ()67p . which can be 
re written as 

{dC,uj) = {C,(ko), 

where dC is the boundary of the p— chain C oriented coherently with C on M. While the 
boundary operator 9 is a global operator, the coboundary operator d is local, and thus more 
suitable for applications. The main property of the exterior differential, 

do d = cF' = =^ d o d = d"^ = 0, (and converse), 

can be easily proved using the Stokes' theorem as 

= {d'^C,oj) = {dC,duj) = {C,d^oj) = 0. 

5.4 De Rham cohomology 

In the Euclidean 3D space we have the following de Rham cochain 

grad curl div 

Using the closure property for the exterior differential in R'^, d o d ^ 
standard identities from vector calculus 

curl • grad = and div • curl = 0. 

As a duality, in we have the following chain complex 

^ Co(M^)^Ci(m2)^C2(m3)^^C3(M^) ^ 0, 
(with the closure property d o d = = 0) which implies the following three boundaries: 

Ci^Co = d{Ci), C2^Ci = d{C2), C3^C2 = d{C3), 



complex 
0. 

= d"^ = 0, we get the 
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where Co G Co is a 0-boundary (or, a point), Ci G Ci is a 1-boundary (or, a line), C2 G C2 is 
a 2~boundary (or, a surface), and C3 G C3 is a 3-boundary (or, a hypersurface). Similarly, 
the de Rham complex implies the following three coboundaries: 

C° = d{C°), ^ C"^ = d{C^), = d{C'^), 

where G Q-^ is 0-form (or, a function), G 0^ is a l~form, G 0^ is a 2-form, and 
G il^ is a 3-form. 

In general, on a smooth nD manifold M we have the following de Rham cochain complex 

m 

^ r?°(M) Jl^(M) ^ r?2(M) f]3(M) ^ ^ ir{M) ^ 0, 

satisfying the closure property on M, do d = d'^ = 0. 

exact forms |?- forms closed forms 




Figure 4: A small portion of the de Rham cochain complex, showing a homomorphism of 
cohomology groups. 

Informally, the de Rham cohomology is the (functional) space of closed differential p— forms 
modulo exact ones on a smooth manifold. 

More precisely, the subspace of all closed p— forms (cocycles) on a smooth manifold M is 
the kernel Ker(d) of the de Rham d— homomorphism (see Figure U]), denoted by ZP(M) C 
QP{M), and the sub-subspace of all exact p— forms (coboundaries) on M is the image Im.{d) 
of the de Rham homomorphism denoted by B'^{M) C ZP{M). The quotient space 

HP (M) = ^'^^^ = Ker {d : nPjM) ^ QP+^M)) 
^ ■ BPM Im(d : QP-^M) ^ nP{M)) ' 
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is called the pth de Rham cohomology group of a manifold M . It is a topological invariant of a 
manifold. Two p— cocycles G r2P(M) are cohomologous, or belong to the same cohomology 
class [a] e HP{M), if they differ by a (p - l)-coboundary a - (5 = dO ^ Vlp-^{M). The 
dimension 6p = dim iJP(M) of the de Rham cohomology group H^j^{M) of the manifold M 
is called the Betti number hp. 

Similarly, the subspace of all p— cycles on a smooth manifold M is the kernel Ker((?) of the 
9— homomorphism, denoted by Zp{M) C Cp{M), and the sub-subspace of all p— boundaries on 
M is the image Im{d) of the 5— homomorphism, denoted by Bp{M) C Cp{M). Two p— cycles 
Ci,C2 G Cp are homologous, if they differ by a (p — 1)— boundary Ci — C2 = OB £ Cp_i(M). 
Then Ci and C2 belong to the same homology class [C] G Hp{M), where Hp{M) is the 
homology group of the manifold M, defined as 

H (M) = ^P(^) = Ker(a:C,(M)^Cp_i(M)) 
>■ Bp{M) lm{d : Cp+i{M) ^ Cp{M)) ' 

where Zp is the vector space of cycles and Bp C Zp is the vector space of boundaries on M. 
The dimension bp = dim Hp{M) of the homology group Hp{M) is, by the de Rham theorem, 
the same Betti number bp. 

If we know the Betti numbers for all (co)homology groups of the manifold M, we can 
calculate the Euler-Poincarc characteristic of M as 

n 

x(M) = Y.{-iYbp. 
p=l 

For example, consider a small portion of the de Rham cochain complex of Figure H] 
spanning a space-time 4-manifold M, 

np-\M) ^ np{M) ^ np+\M) 

As we have seen above, cohomology classifies topological spaces by comparing two subspaces 
of ^l^: (i) the space of p— cocycles, ZP{M) = Kerdp, and (ii) the space of p— coboundaries, 
BP{M) = Imdp-i. Thus, for the cochain complex of any space-time 4-manifold we have, 

^2 = BP{M) C ZP{M), 

that is, every p— coboundary is a p— cocycle. Whether the converse of this statement is true, 
according to Poincare lemma, depends on the particular topology of a space-time 4-manifold. 
If every p— cocycle is a p— coboundary, so that B^ and Z^ are equal, then the cochain complex 
is exact at 0^(M). In topologically interesting regions of a space-time manifold M, exactness 
may fail [37] , and we measure the failure of exactness by taking the pth cohomology group 

HP{M) = ZP{M)/BP{M). 
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